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I.  INTRODUCTION 


This  report  documents  a set  of  analog  operators  written  as  macros 
to  be  used  in  the  Advanced  Continuous  Simulation  Language  (ACSL) . ThfeSfi 
macros  are  compatible  with  the  Electronic  Associates,  Incorporated  (EAI) 
ECSSL  hybrid  compiler. 

The  Systems  Simulation  Directorate  of  MIRADCOM  has  a contract  with 
EAI  to  expand  ECSSL' s capability.  A part  of  the  expansion  is  to  provide 
ECSSL  with  the  ability  to  directly  generate  an  analog  program  design 
from  an  ACSL  digital  computer  program.  That  is,  ECSSL  can  take  an 
ACSL  program  as  its  literal  input  and  generate  an  analog  computer 
representation  of  the  ACSL  program.  To  ensure  compatibility,  ECSSL  must 
be  able  to  generate  analog  representations  of  all  the  pertinent  ACSL 
operators.  ECSSL  also  recognizes  analog  operators  that  do  not  exist  in 
the  standard  versions  of  ACSL.  The  ACSL  macros  discussed  in  this 
document  are  implementations  of  these  special,  analog-oriented  operators. 


[ Many  of  the  algorithms  for  these  operators  are  in  the  original  form 

used  by  EAI  [1];  however,  the  counters,  limited  mode  controlled  integrator, 
and  resolver  macros  have  major  changes  while  still  other  operators  have 
minor  changes.  Ed  Mitchell  and  Joe  Gauthier  of  Mitchell  and  Gauthier 
Associates,  the  authors  of  ACSL,  have  been  of  major  assistance  in  both 
developing  the  algorithms  and  interfacing  the  algorithms  with  ACSL. 

1 Of  course  the  responsibility  for  the  proper  operation  of  these  macros 

lies  with  the  writer  of  this  report. 


The  following  paragraphs  contain  a functional  description  of  each 
operator,  the  reasoning  behind  the  algorithm  for  an  operator  where  an 
explanation  is  believed  to  be  appropriate,  flow  diagrams  where  needed 
and  a description  of  the  verification  tests  for  each  algorithm.  There 
is  also  a general  discussion  section  describing  certain  aspects  of  macro 
design,  implementation  and  usage. 


II.  GENERAL  DISCUSSION 


There  are  a number  of  special  programming  rules  that  should  be 
followed  when  using  or  programming  these  macros.  (Macro  programming  in 
general  is  discussed  in  Reference  [21.)  Comments  applicable  to  the 
macros  in  general  are  mentioned  here.  Restrictions  peculiar  to  a 
particular  macro,  such  as  regions  of  validity,  integration  step  size, 
etc.,  are  described  in  the  section  for  that  macro.  Special  ACSL  internal 
functions  and  variables  needed  in  the  macro  code  are  elaborated  upon  in 
this  section.  System  control  statements  needed  to  run  an  ACSL  program 
on  the  CDC  6600  and  CYBER  74  are  shown. 


A.  Macro  Implementation  and  Usage 


In  the  following  sections  where  a listing  of  macro  coding  on 
cards  is  shown,  it  becomes  apparent  that  the  code  has  been  pushed  to  the 
left  side  of  the  cards,  that  is,  the  coding  starts  in  Column  1 and  most 
of  the  code  contains  no  blanks.  This  was  done  to  minimize  the  amount  of 
computer  memory  needed  to  store  the  macros  during  program  compilation 
and  execution.  For  the  same  reason,  there  are  no  comment  cards.  The 
memory  allocation  of  the  user's  ACSL  program  code  is  automatically 
minimized  by  the  ACSL  compiler. 

An  input  variable  in  a macro  argument  list  can  take  the  form  of  a 
constant,  a logic  expression  or  arithmetic  expression,  as  well  as,  a 
variable  name.  ACSL  takes  the  macro  input  list  as  invoked  by  the  user 
and  inserts  this  in  replacement  of  the  original  associated  variable 
names  in  the  macro  program  code.  This  occurs  when  the  user's  program 
code  is  being  compiled  into  FORTRAN  code.  For  example,  suppose  an  ACSL 
program  is  as  follows: 

PROGRAM  ACSL  EXAMPLE 


BADCTR  (LYO,  LY1 , LY2,  LY3,  LYCO 


I+J-l,  I.XUP,  LXCI ) 


END  $ ’OF  PROGRAM' 


where  the  macro  called  is 

MACRO  BADCTR  (LYO, LY1 ,LY2,LY3, LYCO, NIC, LXUP, LXCI) 


K = (NIC) *10 


MACRO  END 

The  user  has  inserted  an  arithmetic  expression  for  NIC.  The  FORI  RAN 
code  generated  by  the  ACSL  compiler  will  replace  the  macro  code 

K * (NIC) *10 
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by 


m 


ZZ...  = (K+J-l) *10 

where  ZZ...  is  the  global  name  assigned  to  the  local  macro  variable  K 
(A  local  variable  is  not  an  input  or  output  of  the  macro.).  It  should 
be  noted  that  all  local  variables  must  appear  in  a MACRO  REDEFINE 
statement . 

Creating  macros  out  of  the  algorithms  required  the  use  of  the  ACSL 
system  variables  ZZICFL  and  ZZRNFL  and  the  ACSLLIB  function 

ZZFST  (DUMMY) . 

ZZICFL  is  a logical  variable  that  is  set  to  "TRUE"  when  the  ACSL 
Executive  is  in  the  INITIAL  section  of  an  ACSL  program,  otherwise  it  is 
set  to  "FALSE".  ZZICFL  is  used  to  determine  the  appropriate  time  to  set 
initial  conditions.  ZZRNFL  is  a logical  variable  that  is  set  to  "FALSE" 
except  between  the  time  the  run  time  command  REINIT  is  invoked  and  the 
time  the  ACSL  Executive  finishes  execution  of  the  INITIAL  section. 

ZZRNFL  is  used  to  initialize  the  macro  variables  to  the  values  existing 
at  the  end  of  the  last  run.  ZZFST  is  an  ACSLLIB  FORTRAN  real  function 
whose  argument  is  a dummy  variable.  It  essentially  has  a value  of  1.0 
during  the  first  intermediate  integration  step  of  a multistep  integration 
algorithm  and  has  approximately  a value  of  0.0  otherwise.  All  macro 
algorithms  that  are  executed  only  once  per  major  integration  step  do  so 
when  ZZFST  has  a value  of  1.0.  For  example,  the  logic  macros  are 
executed  only  once  per  major  integration  step.  This  corresponds  to  the 
analog  logic  operations  changing  states  once  per  analog  clock  cycle. 
Because  ZZFST  is  a FORTRAN  function,  it  can  always  be  used  despite 
internal  changes  to  ACSL.  However,  the  two  logical  variable  names 
ZZICFL  and  ZZRNFL  are  subject  to  ACSL  compiler  changes.  If  changed, 
they  would  cause  erroneous  results.  It  should  be  noted  that  these 
specific  macros  are  the  only  ACSL  programming  features  available  to  the 
user  that  could  be  affected  in  this  manner.  It  has  also  been  assumed 
(and  is  true  in  the  present  version  of  ACSL)  that  when  ZZRNFL  is  "TRUE" 
during  a call  from  the  INITIAL  section  to  the  DERIVATIVE  section  of  the 
program  ZZICFL  is  also  "TRUE". 

The  FORTRAN  logical  operator  AND  can  be  used  with  integer  variables 
and  integer  constants  under  the  CDC  SCOPE  and  NOS/BE  systems.  This  is 
not  the  case  with  many  other  computers.  But  ACSLLIB  has  its  own  FORTRAN 
integer  function 

AND  (A, B) 

which  contains  an  installation  dependent  algorithm  that  can  perform  an 
AND  with  integers.  Hence,  ACSL  macros  in  Sections  IV  and  V that  use 
AND  with  integers  should  function  on  any  computer  system  using  ACSL. 
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Output  argument  variables  that  are  logical  are  so  declared  inside 
the  macro.  However,  if  the  user  calls  the  macro  in  the  form 

LY=MODE(LX) 
rather  than 

MODE (LY=LX) , 

LY  must  also  be  declared  logical  in  the  user's  program.  Good  programming 
practice  dictates  that  all  known  logical  variables  in  the  user's  program 
be  so  declared  by  the  user.  By  declaring  logical  output  variables  inside 
a macro,  a programmer  ensures  that  if  one  macro  calls  another  macro  that 
has  a logical  output  argument,  which  is  essentially  unseen  by  the  user, 
then  the  argument  will  be  properly  declared. 

When  the  ACSI.  run  time  command  REINTT  is  invoked,  no  changes  are 
made  to  any  of  the  outputs  of  a macro.  The  initial  condition  code  is 
skipped  and  the  correct  transition  from  the  last  simulation  run  to  the 
next  run  results.  However,  these  macros  will  not  be  properly  initialized 
for  a third  run  where  a REINIT  command  does  not  precede  a START  command. 
For  example,  a run  time  command  sequence  will  run  correctly,  as  follows 

START. . .REINIT. . .START. . .REINIT. . .START 
but  the  third  run  in  the  following  run  time  command  sequence  will  not, 

START. . .REINIT. . .START. . .START 

The  difficulty  is  that  special  consideration  must  be  given  to  analog 
hardware  since  it  does  not  run  in  a REINIT  mode.  Digitally  it  is 
possible  to  make  the  analog  macros  initialize  correctly  under  REINIT 
conditions,  but  not  without  proper  consideration  by  the  user  of  burden- 
some details  the  macros  are  specifically  designed  to  relieve.  In  sum- 
mary, the  outputs  of  the  analog  macros  will  transition  properly  for  the 
situation  where  a REINIT  command  is  used  between  two  simulation  runs, 
but  the  initial  conditions  on  the  macros  will  not  change  to  the  proper 
new  values  to  allow  follow-on  runs  that  are  not  preceded  by  a REINIT 
statement.  The  one  exception  to  this  is  the  I.MMINT  macro  in  Section  XII. 
I.MMINT  functions  properly  under  any  run  time  command  sequence  containing 
a REINIT. 

Analog  logic  usually  contains  algebraic  loops.  Consequently  the 
resulting  ACSL  code  for  such  analog  logic  will  have  to  be  embedded  in 
PROCEDURAL  statements.  It  is  not  always  a straightforward  task  to  get 
the  ACSL  compiler  to  accept  such  code.  An  example  is  shown  in  Section  IV 
where  two  BCDCTR  statements  are  embedded  in  a PROCEDURAL.  Embedding 
just  one  of  them  in  a PROCEDURAL  was  not  acceptable  to  the  ACSL  compiler. 


B.  Interactive  anil  Batch  .lob  Control  Statements 


Figure  1 contains  the  NOS/BE  interactive  terminal  commands 
for  the  CYBER  74  system  which  will  compile  and  execute  an  ACSL  program. 
Figure  2 contains  the  NOS/BE  interactive  terminal  commands  for  the 
CYBER  74  system  which  will  compile  an  ACSI.  program  and  generate  an 
absolute  file  which  in  turn  will  produce  hard  copy  plots  on  the 
Tektronix  4014.  Figure  3 contains  the  SCOPE  job  control  cards  to  run 
batch  jobs  on  the  CDC  6600. 

'R, INPUT'  in  Figures  I and  2 is  the  interactive  command  to  read  the 
user's  ACSI,  program  from  the  card  reader  and  name  the  resulting  local 
file  'INPUT'. 


III.  AMODPI 

This  macro  is  used  by  the  analog  resolver  macros  to  output  an 
angle  which  lies  in  the  interval  |-  , + i | given  any  input  angle.  It  is 
also  used  to  provide  a logical  variable  that  can  generate  an  angle  as  it 
would  be  output  from  a coni  imious  resolver.  The  standard  torm  is 

AMODPI  ( Tll,l,  Y IMS  ) 

where  TH  is  an  angle  (radians)  between  |-i,+n|,  I.Y  is  a logical  variable 
and  TUS  is  the  input  angle  (radians).  Figure  4 illustrates  the  above 
for  the  case  where  TUS  is  a ramp  in  I imo.  TDK  is  the  output  as  would 
be  seen  from  a cent  imious  resolver.  From  Figure  4 one  can  see  that  THR 
can  easily  be  generated  from  the  AMODPI  outputs  hv  flipping  the  sign  on 
TH  when  I.Y  is  FAI.SE.  To  derive  the  algorithm  used  for  AMODPI  the  fol- 
lowing, should  be  noted  from  Figure  4: 

a)  For  TUS  0,  I.Y  FAI.SE,  ■ TIIS  In,  5tt<TIIS<7m,  etc. 

TRUE,  0 TIIS_  n , In  TIIS  5n  , etc . 

b)  When  I.Y  is  TRUE,  (TIIS-  I'll)  / 2 1 is  an  even  integer.  When  I.Y  is 
FAI.SE,  ( TIIS— Til ) / 2 i is  an  odd  integer. 

c)  Therefore  when  III  is  lotnul,  I.Y  can  be  determined  from  b)  . 

For  TUS- 0 the  criterion  for  assigning  a value  oi  TRUE  or  FAI.SE  to  I.Y 
still  holds. 


Figure  S illustrates  the  angles  used  in  the  I ol lowing  equation  for 


TH 


III  = < 


2 « * I NT  (a/ii) 


* .E.V  ‘ 


where 


a = THS  (modulo  2tt) 

and  INT  determines  the  integer  portion  of  its  argument.  This  can  also 
be  written  as 

TH  = 6 - tt  * INT(a/n) 


where 


6 = THS  (modulo  tt)  . 

Using  Figure  5 and  knowing  that 

a = THS(modulo  2tt)  it  can  be  determined  that 
for  THS  > 0, 

INT(THS/tt)  = odd  integer  in  III  and  IV 

even  integer  in  I and  II  and 

for  THS  < 0, 

INT(THS/tt)  = even  integer  in  III  and  IV 
odd  integer  in  T and  II. 

Therefore 

TH  * THS  (modulo  tt)  - ti*[  INT(THS/Tt)  (modulo  2)]. 
If  N is  defined  as  follows, 

N = INT  (THS/ tt) 


then 


TH  = THS  - N*tt  - tt  (N-2*INT(N/2) ) 

= THS  - 2tt  * (N-  INT (N/2)  ) . 

It  has  been  stated  from  observations  of  Figure  4 that  LY  should  be 
TRUE  or  FALSE  depending  on  whether  (THS-TH)/2ti  is  even  or  odd. 

From  above 

(THS-TH)  / 2 tt  = N - INT  (N/2)  . 

Therefore  LY  can  be  found  by  the  following 

[N-INT(N/2)]  (modulo  2)  = 0,  LY  = TRUE 

1,  LY  = FALSE 


Figure  6 is  a listing  of  the  AMODPI  macro.  Note  that  the  intrinsic 
functions  MOD(...)  and  INT(...)  must  be  declared  integer.  No  testing 
was  done  of  this  macro  alone,  but  it  was  checked  out  through  the  tests 
for  the  resolvers. 


IV.  BCDCTR 

The  macro  models  a binary  coded  decimal  counter  with  a A-bit 
register.  The  standard  form  is 

BCDCTR (I.YO , I.Y1 , I.Y2,  LY3,  I.YCO  = LXl,  I.X2,  I.X3,  l.XUP,  LXCI,  RST)  . 

LYO  through  LY3  are  output  logic  variables  whose  values  (TRUE  or  FALSE) 
correspond  to  the  states  (high  or  low)  of  the  four  flip-flops  in  the 
register.  LYO,  LY1,  I.Y2,  and  LY3  represent  1,  2,  A,  and  8,  respectively. 
I.XO,  LX1,  LX2 , LX  1 are  logical  variables  or  expressions  whose  values  cor- 
respond to  tiie  initial  states  of  the  four  flip-flops  in  the  register  or 
to  the  states  into  which  the  flip-flops  change  when  the  logical  variable 
or  logical  expression  RST  becomes  TRUE.  The  latter  is  essentially  the 
reset  operation  capability  many  analog  machines  have.  l.XUP  and  LXCI 
are  logical  input  variables  or  expressions.  If  l.XUP  is  TRUE  the  regis- 
ter adds  one  to  the  count  for  every  major  integration  step  during  wtiich 
LXCI  is  TRUE.  If  l.XUP  is  FALSE  the  register  subtracts  one  from  the 
count  on  everv  major  integration  step  during  which  LXCI  is  TRUE.  The 
logical  output  variable  I.YCO  is  the  carry-over  output.  It  will  be  TRUE 
for  one  major  integration  step  whenever  the  register  is  counting  up  and 
contains  nine  (bit  pattern  1001)  or  when  the  register  is  counting  down 
and  contains  zero  (bit  pattern  0000).  Thus  registers  can  be  cascaded 
by  connecting  the  CO  of  one  register  to  the  Cl  of  the  next. 

Figure  7 contains  a I low  diagram  ol  BCDCTR.  Figure  8 is  an  ACSL 

listing  of  BCDCTR  where  the  AND  ( ) is  an  ACSLLIB  FORTRAN 

function  and  must  be  declared  integer.  Referencing  the  flow  diagram, 
if  7.7.FST(K)  is  less  than  0.5  or  Z7RNFI.  is  TRUE,  then  the  rest  of  the 
macro  code  is  bypassed.  Hence,  this  macro  is  performed  only  once  per 
major  integration  step  and  when  the  REIN1T  run  time  command  is  invoked, 
no  values  are  changed.  If  77ICKL  is  TRUK,  then  K,  the  present  integer 
number  representation  in  the  register,  is  set  to  the  equivalent  value 
represented  by  1.X0  through  LX  1 and  I.III,  the  local  variable  for  I.YCO,  is 
initialized  to  FALSE.  The  macro  then  jumps  to  a part  of  the  code  that 
is  explained  below.  II  local  logical  variable  NABI.E  is  FALSE,  the  macro 
jumps  to  the  last  line  of  code,  which  sets  NABI.E  to  the  value  in  LXCI. 
This  causes  NABI.E  to  become  TRUE  one  major  integration  step  after  LXCI 
first  becomes  TRUE.  The  one  integration  step  lag  is  purposefully  added 
to  ensure  proper  counting  when  registers  are  cascaded  together.  When 
RST  is  TRUE  the  program  branches  as  when  ZZICFL  is  TRUE.  However, 
setting  LHI  to  FALSE  when  RST  is  TRUE  does  something  special.  It  ensures 


r - 

that  LYCO  is  set  TRUE  (instead  of  LYO  through  LY3  and  K being  reset 
again  after  RST  has  already  caused  them  to  be  reset)  if  LXUP  is  TRUE 
and  K is  9 or  if  LXUP  is  FALSE  and  K is  0.  If  LXUP  is  TRUE,  K is 
incremented  by  one;  if  LXUP  is  FALSE,  K is  decremented  by  one.  KP  is 
the  previous  decimal  representation  in  the  register.  LYCO  is  set  FALSE 
to  keep  LYCO  from  being  TRUE  for  more  than  one  major  integration  step 
(clock  pulse)  while  LXCI  is  TRUE.  If  K is  equal  to  0 and  LXUP  is  FALSE, 
or  if  K is  equal  to  9 and  LXUP  is  TRUE;  then  LHI  is  checked.  If  LHI  is 
FALSE  LYCO  is  set  TRUE;  if  LHI  is  TRUE  LYO  through  LY3  are  set  to  rep- 
resent 0 or  9 and  K is  set  to  0 or  9,  depending  on  whether  LXUP  is 
FALSE  or  TRUE,  respectively.  LHI  is  set  equal  to  LYCO  later  in  the  code, 
thus  causing  a lag  of  one  major  integration  step  between  LYCO  and  LY. 

This  ensures  that  LYO  through  LY3  and  K are  reset  the  integration  step 

after  K is  equal  to  0 or  9.  If  K is  neither  0 nor  9 then  LYO  through 
LY3  are  calculated  by  AND'ing  K with  integers  1,  2,  4,  and  8,  respec- 
tively. KP  is  set  to  K,  LHI  is  set  to  LYCO,  and  NABI.E  is  set  to  LXCI. 

Note  that  LYCO  cannot  be  substituted  into  the  argument  list  for  RST 
when  invoking  the  BCDCTR  or  BINCTR  macro.  However,  RST  can  be  set  equal 
to  LYCO  in  a logical  expression. 


Figure  9 is  a listing  of  the  ACSL  test  program  for  BCDCTR.  Two 
BCD  registers  are  cascaded  together.  The  reset  feature  of  the  macro  is 
also  exhibited.  The  test  ensures  that  the  macro  counts  up  and  down 
properly  with  LYCO  (actually  I.XH  and  LYH)  switching  at  the  appropriate 
time.  Also  checked  is  the  functioning  of  the  macro  when  the  REINIT  run 
time  command  is  invoked.  The  BCDCTR  macro  is  called  twice  in  the 
DERIVATIVE  section  and  represents  two  BCD  registers  cascaded  together. 
The  logical  variables  LXO  through  LX3  represent  a number  in  the  ones 
unit  place,  and  variables  LYO  through  LY3  represent  a number  in  the  tens 
unit  place.  The  constants  AO  througli  A3  and  BO  through  B3  defined  in 
the  INITIAL  section  are  the  initial  conditions  and  reset  values  for 
LXO  through  I.X3  and  LYO  through  LY3,  respectively.  AO  through  A3  repre- 
sent 8 and  BO  through  B3  represent  80.  LXUP  and  LXCI  are  functions  of 
time  as  defined  in  the  DERIVATIVE  section  via  the  I.SW  and  CLOCK 
(Section  VI)  macros.  LXH,  the  carry-over  output  bit  for  the  ones  place 
number,  is  AND'ed  with  LXCI  to  form  the  carry-over  enablement  input  to 
the  tens  place  number.  LYH  is  the  carry-over  output  for  the  tens  place 
number.  LHICH  is  the  logical  input  to  both  registers.  In  its  equation 
in  the  DERIVATIVE  section  LHICH  becomes  TRUE  when  either  the  two 
counters  contain  99  and  are  counting  up  or  the  two  counters  show  00  and 
are  counting  down.  Whenever  LHICH  becomes  TRUE  the  cascaded  counters 
are  reset  to  88.  Note  that  the  cascaded  counters  and  the  LHICH  expres- 
sion are  embedded  in  a PROCEDURAL  (where  LHICH  is  the  input  and  LYH  is 
the  output)  due  to  an  algebraic  loop.  LHICH  is  calculated  outside 
the  BCDCTR  argument  lists  because  it  contains  LXH.  In  the  DYNAMIC 
section,  NUMBER,  the  numerical  value  for  the  cascaded  counters,  is 
calculated  for  display  purposes. 


10 


. i 


In  the  run  time  commands  the  variables  to  be  printed  are  listed  in 
a PREPAR  statement.  Dump  is  set  to  TRUE  so  a debug  printing  is  made 
from  the  TERMINAL  section.  TSTP,  the  simulation  stopping  time,  is  set 
to  40  sec  and  NDBIK'  is  set  to  5 so  that  most  of  the  variables  for  the 
first  five  integrations,  counting  minor  integration  steps,  are  printed. 
PROCED  'GO'  when  invoked  will  start  the  simulation  and  print  the  data 
according  to  the  PREPAR  list.  The  RE1N1T  command  is  invoked,  TSTP  is 
set  to  240  sec  and  PROCED  'GO'  is  again  invoked.  Figure  10  shows  the 
printed  results.  The  variable  NUMBER  lias  repeated  values  due  to  the 
CLOCK  input  in  the  equation  for  LXC1 . 


V.  BINCTR 

,(  The  macro  models  a binary  counter  with  a 4-bit  register.  The 

standard  form  is 

BINCTROYO,  1.Y1.  1.Y2,  LY3,  LYCO  = LXO,  LX1 , LX2 , LX3,  LXUP, 

LXC1,  RST) . 

LYO  through  LY3  are  output  logical  variables  whose  values  (TRUE  or 
FALSE)  correspond  to  the  states  (high  or  low)  of  the  four  flip-flops  in 
the  register.  LYO,  LY 1 , I.Y2,  and  l.Y  ) represent  1,  2,  4,  and  8 respec- 
tively. LXO,  LX1 , LX2,  and  1.X3  are  logical  variables  or  expressions 
whose  values  correspond  to  the  initial  states  of  the  four  flip-flops  in 
the  register  or  the  states  t lie  flip-flops  change  to  when  the  logical 
variable  or  expression  RST  becomes  TRUE.  The  latter  is  essentially  the 
reset  operation  capability  many  analog  machines  have.  LXUP  and  LXCI 
are  logical  input  variables  or  expressions.  If  LXUP  is  TRUE  the  register 
.adds  one  to  the  count  every  major  integration  step  during  which  LXCI  is 
TRUE.  If  LXUP  is  FALSE  the  register  subtracts  one  from  the  count  on 
every  major  integration  step  during  which  LXCI  is  TRUE.  The  logical 
output  variable  LYCO  is  the  carry-over  output.  It  will  become  TRUE  for 
one  major  integration  step  whenever  the  register  is  counting  up  and 
contains  fifteen  (bit  pattern  1111)  or  when  the  register  is  counting 
down  and  contains  zero  (bit  pattern  0000).  Thus  registers  can  ' • 
cascaded  by  connecting  the  CO  of  one  register  to  the  Cl  of  t lie  next . 

Figure  11  contains  a flow  diagram  of  BINCTR  and  Figure  12  is  an 
ACSI.  listing  of  BINCTR.  Section  IV,  the  BCDCTR  description,  has  a 
narrative  on  the  flow  diagram.  The  two  diagrams  are  essentially  the 
same . 

Figure  13  is  the  BINCTR  test  program  listing  - it  checks  the 
operation  of  just  one  register.  The  test  primarily  ensures  the  macro 
counts  up  and  down  properly  with  LYCO  switching  at  the  appropriate  time. 
The  reset  capability  is  also  tested. 
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The  BINCTR  macro  is  called  in  the  DERIVATIVE  section.  The  con- 
stants BO  through  B3  defined  in  the  INITIAL  section  are  the  initial 
conditions  and  reset  values  for  LYO  through  LY3,  respectively.  BO 
through  B3  represent  8.  LXUP,  I.XCI  and  LRST  are  functions  of  time  as 
defined  in  the  DERIVATIVE  section  via  the  I.SW  macro.  In  the  DYNAMIC 
section,  NUMBER,  the  numerical  value  in  the  register,  is  calculated  for 
display  purposes. 

In  the  run  time  commands  the  variables  to  be  printed  are  listed  in 
a PREPAR  statement.  TSTP,  the  simulation  stopping  time,  is  preset  in 
the  INITIAL  section  to  38  sec.  The  first  ten  integration  steps  are 
printed  for  most  of  the  program  variables  by  setting  NDBUG  = 10. 

PROCF.D  'GO'  when  invoked  will  start  the  simulation  and  print  the  data 
according  to  the  PREPAR  list.  The  REINIT  command  is  invoked,  TSTP  is 
set  to  45  sec  and  'GO'  is  again  invoked.  Figure  14  shows  the  printed 
data  from  the  simulation  run.  Note  that  for  T = 3.0  sec  LXCI  and  LXUP 
go  TRUE  and  FALSE,  respectively,  in  the  middle  step  of  the  RK2  inte- 
gration algorithm  and  are  printed  at  the  end  of  an  integration  step. 
However  LYO  through  I.Y3  and  l,YC0  change  values  only  during  the  first 
step  of  a RK2  algorithm  so  their  printed  values  will  show  a change  two 
integration  steps  after  the  one  where  I.XCI  and  LXUP  has  a printed 
value  change. 


VI.  CLOCK 

The  macro  models  a clock  pulse  train.  The  standard  form  is 
LY  = CLOCK  (R). 

LY  is  a logical  variable  and  R is  usually  a real  constant  defining  the 
time  step  between  pulses.  I.Y  is  defined  as  follows: 

I.Y  = .F.  for  the  independent  variable  (usually  time) 

i nR,  n = 0,1,2... 

LY  = .T.  for  the  independent  variable  = nR 

The  macro  is  listed  in  Figure  15.  The  clock  pulses  are  generated  by 
calling  the  differentiator  macro  (Section  X)  with  the  LPULSE  macro 
(Section  XIII)  as  the  input  argument. 

Figure  16  is  the  ACSL  test  program  listing  for  the  CLOCK,  LPULSE, 
MODE,  SWITCH  and  I.STF.P  macros.  CLOCK  is  invoked  in  the  DERIVATIVE 
section  with  its  output  variable  being  LY 1 and  its  input  being  PERIOD, 
which  is  set  to  2.0  sec  by  a CONSTANT  statement  in  the  INITIAL  section 


In  the  run  time  commands  the  variables  to  be  printed  are  listed  in 
a PREPAR  statement.  Dump  is  set  to  TRUE  so  a debug  printing  is  made 
from  the  TERMINAL  section.  TSTP,  the  simulation  stopping  time,  is  set 
to  AO  sec  and  NDBUG  is  set  to  5 so  that  most  of  the  variables  for  the 
first  five  integrations,  counting  minor  integration  steps,  are  printed. 
PROCEI)  'GO'  when  invoked  will  start  the  simulation  and  print  the  data 
according  to  the  PREPAR  list.  The  RE  IN  IT  command  is  invoked,  TSTP  is 
set  to  2A0  sec  and  PROCED  'CO'  is  again  invoked.  Figure  10  shows  the 
printed  results.  The  variable  NUMBER  has  repeated  values  due  to  the 
('LOCK  input  in  the  equation  for  LXCI. 

V.  BINCTR 

The  macro  models  a binary  counter  witli  a 4-bit  register.  The 
standard  form  is 

BINCTRfLYO,  LY1,  LY2 , I.Y3,  LYCO  = LXO,  LX1  , LX2 , 1.X3 , LXUP, 

LXCI,  RST) . 

LYO  through  LY3  are  output  logical  variables  whose  values  (TRUE  or 
FALSE)  correspond  to  the  states  (high  or  low)  ot  the  four  flip-flops  in 
the  register.  LYO,  I.Y1,  LY2,  and  l.Y  1 represent  1,  2,  4,  and  8 respec- 
tively. LXO,  LX1,  I.X2,  ami  LX 3 are  logical  variables  or  expressions 
whose  values  correspond  to  the  initial  states  of  the  four  flip-flops  in 
the  register  or  the  states  the  I lip-1  lops  change  to  when  the  logical 
variable  or  expression  RST  becomes  TRUE.  The  latter  is  essentially  the 
reset  operation  capability  many  analog  machines  have.  LXUP  and  LXCI 
are  logical  input  variables  or  expressions.  II  LXUP  is  TRUE  the  register 
adds  on  to  the  count  every  major  integration  step  during  which  LXCI  is 
TRUE.  If  LXUP  is  FALSE  the  register  subtracts  one  from  the  count  on 
every  major  integration  step  during  which  LXCI  is  TRUE.  The  logical 
output  variable  LYCO  is  the  carry-over  output.  It  will  become  TRUE  for 
one  major  integration  step  whenever  the  register  is  counting  up  and 
contains  fifteen  (bit  pattern  1111)  or  when  the  register  is  counting 
down  and  contains  zero  (bit  pattern  0000).  Thus  registers  can  be 
cascaded  by  connecting  the  CO  ol  one  register  to  the  Cl  of  the  next. 

Figure  11  contains  a I low  diagram  ol  BINCTR  and  Figure  12  is  an 
ACSL  listing  of  BINCTR.  Section  IV,  the  BCDCTR  description,  has  a 
narrative  on  the  flow  diagram.  The  two  diagrams  are  essentially  the 
same . 


Figure  13  is  the  BINCTR  test  program  listing  - 
operation  of  just  one  register.  The  test  primarily 
counts  up  and  down  properlv  with  LYCO  switching  at 
The  reset  capability  is  also  tested. 


it  checks  the 
ensures  the  ma.ro 
the  appropriate  time. 


The  BINCTR  macro  is  called  in  the  DERIVATIVE  section.  The  con- 
stants BO  through  B3  defined  in  the  INITIAL  section  are  the  initial 
conditions  and  reset  values  for  LYO  through  LY3,  respectively.  BO 
through  B3  represent  8.  LXUP,  LXCI  and  LRST  are  functions  of  time  as 
defined  in  the  DERIVATIVE  section  via  the  l„SW  macro.  In  the  DYNAMIC 
section,  NUMBER,  the  numerical  value  in  the  register,  is  calculated  for 
display  purposes. 

In  the  run  time  commands  the  variables  to  be  printed  are  listed  in 
a PREPAR  statement.  TSTP,  the  simulation  stopping  time,  is  preset  in 
the  INITIAL  section  to  38  sec.  The  first  ten  integration  steps  are 
printed  for  most  of  the  program  variables  by  setting  NDBUG  = 10. 

PROCED  'GO'  when  invoked  will  start  the  simulation  and  print  the  data 
according  to  the  PREPAR  list.  The  REINIT  command  is  invoked,  TSTP  is 
set  to  45  sec  and  'GO'  is  again  invoked.  Figure  14  shows  the  printed 
data  from  the  simulation  run.  Note  that  for  T = 3.0  sec  LXCI  and  LXUP 
go  TRUE  and  FALSE,  respectively,  in  the  middle  step  of  the  RK2  inte- 
gration algorithm  and  are  printed  at  the  end  of  an  integration  step. 
However  LYO  through  LY3  and  LYCO  change  values  only  during  the  first 
step  of  a RK2  algorithm  so  their  printed  values  will  show  a change  two 
integration  steps  after  the  one  where  l.XCI  and  LXUP  has  a printed 
value  change. 


VI.  CLOCK 

The  macro  models  a clock  pulse  train.  The  standard  form  is 
LY  = CLOCK  (R). 

LY  is  a logical  variable  and  R is  usually  a real  constant  defining  the 
time  step  between  pulses.  LY  is  defined  as  follows: 

I.Y  = .F.  for  the  independent  variable  (usually  time) 

i nR,  n = 0,1,2... 

LY  = .T.  for  the  independent  variable  = nR 

The  macro  is  listed  in  Figure  15.  The  clock  pulses  are  generated  by 
calling  the  differentiator  macro  (Section  X)  with  the  l.PULSE  macro 
(Section  XIII)  as  the  input  argument. 

Figure  16  is  the  ACSL  test  program  listing  for  the  CLOCK,  LPULSE, 
MODE,  SWITCH  and  LSTF.P  macros.  CLOCK  is  invoked  in  the  DERIVATIVE 
section  with  its  output  variable  being  LY1  and  its  input  being  PERIOD, 
which  is  set  to  2.0  sec  by  a CONSTANT  statement  in  the  INITIAL  section 
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In  the  run  time  commands  the  data  to  be  printed  are  listed  in  a 
PREPAR  statement.  NDBUG  is  set  to  8.  This  causes  the  printing  of 
most  of  the  simulation  variables  for  the  first  eight  integration  steps 
including  the  minor  integration  steps  tor  a multistep  integration 
algorithm.  TSTP,  the  simulation  stopping  time,  is  set  to  2.0  sec.  The 
START  command  causes  the  actual  running  of  the  simulation.  In  the 
DYNAMIC  section  of  the  program  the  simulation  is  stopped  by  the  TERMT 
macro  when  the  independent  variable,  time,  exceeds  TSTP.  The  PRINT 
command  causes  the  printing  of  the  data  in  the  PREPAR  statement  every 
communication  interval,  C1NT.  TINT  is  defined  to  be  1 sec  in  the 
FNITIAl.  section.  The  RE1N1T  command  is  invoked  and  TSTP  is  reset  to 
8 sec.  The  simulation  is  run  again  and  the  resulting  data  printed. 
Figure  17  contains  the  data  t tom  t lie  two  simulation  runs  as  printed  on 
a line  printer.  The  data  are  for  the  1, PULSE,  MODE,  SWITCH  and  ESTEP 
macros,  as  well  as,  for  the  CLOCK  macro. 


VII  CPTR 

The  macro  models  a continuous  polar  to  rectangular  resolver. 
The  standard  form  is 

CPTR  (X,  Y,  Til,  1 Y R,  THPOT,  I'll  10) 


where 

TH  = f[  THPOT  dt  + I'll  I C (as  would  he  output  from  a resolver) 

( > 

X = R*C0S(TH) 

Y = R*S  IN(Tll) 

I.Y  = .TRUE.  , -n  TH  n , (4n-l)n  I'll  (4nfl)it 

( t-4n ) ti  i'll  (-  l-4n) « 

= .FALSE.,  (An-l)n-TIT  (4n-l)ii 
( )-4n ) n 111(1  -An ) « 

n = a positive  integer. 

All  angles  and  angle  derivatives  are  in  radians  and  radians  per  second. 

Figure  18  is  a I low  diagram  of  CPTR  and  Figure  ID  is  a CPTR  listing. 
CPTR  is  calculated  every  intermediate,  as  well  as,  every  major  integra- 
tion step  so  there  is  no  need  to  check  X'/.FST . TII  as  output  from  a 
resolver  must  lie  between  the  limits  ol  1 n no  matter  how  large  the 
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actual  angle  THS  becomes.  However  the  initial  condition  THIC  should  be 
able  to  be  input  without  the  +ti  restriction.  But  to  obtain  correct 
values  for  LY,  THS  must  be  integrated  with  an  initial  condition,  THICP, 
that  lies  within  ±n.  Hence  when  ZZICFL  is  TRUE  macro  AMODPI  reduces 
THIC  to  a trigonometrically  equivalent  angle  within  in  and  returns  this 
value  in  THICP.  This  is  exactly  what  an  analog  programmer  must  do  to 
set  the  initial  condition  on  a resolver.  When  ZZRNFI.  is  TRUE,  THICP  is 
set  to  the  value  of  THS  at  the  end  of  the  last  run.  THDOT  is  integrated 
to  obtain  the  actual  angle  THS  (except  for  the  difference  between  THIC 
and  THICP).  It  is  to  be  noted  that  in  Figure  19  THICP  is  enclosed  in 
parentheses  in  the  statement 

THS  = INTEC  (THDOT, (THICP)) 

to  let  the  ACSL  Executive  know  the  initial  condition  for  this  state 
variable  must  be  calculated.  The  ACSL  PTR  macro  changes  the  polar 
coordinates  R and  THS  to  the  rectangular  coordinates  X and  Y.  AMODPI 
converts  THS  to  THP,  its  equivalent  value  within  lit.  Using  LY  and 
THP,  TH  can  be  derived  from  the  ACSL  RSW  operator.  In  Section  III  there 
is  a discussion  of  the  process  for  obtaining  TH. 

Figure  20  is  the  CPTR  test  program  list  ing.  CPTR  is  invoked  in 
the  DERIVATIVE  section.  R is  set  to  1.0  in  the  INITIAL  section  by  a 
CONSTANT  statement.  THDOT  and  THIC  are  set  to  1.745  rad/sec  and  0.1745 
rad,  respectively.  LYG  is  a function  of  LY  that  takes  on  the  value  of 
1.0  or  0.0  when  I,Y  is  TRUE  or  FALSE,  respectively.  In  the  DYNAMIC 
section  THC,  THDOTG  and  THICC  arc  defined  as  being  TH , THDOT  and  THIC 
expressed  in  degrees,  where  RADEG  is  the  conversion  factor  from  radians 
to  degrees  as  specified  in  the  INITIAL  section. 

Figure  21  contains  the  run  time  commands  entered  on  the  Tektronix 
4014  to  obtain  plots  and  deferred  printer  output.  The  ACSL  system 
symbols  PRN  and  CMD  are  set  to  9 and  D1S,  respectively.  The  former 
generates  a separate  file  named  'PRINT'  that  contains  all  the  run  time 
commands  entered  on  the  Tektronix  and  all  the  data  generated  by  the 
PRINT  and  DISPLY  run  time  commands.  The  latter  declares  that  the 
logical  unit  on  which  display  data  are  written  is  the  same  logical  unit 
from  which  the  run  time  commands  are  read.  This  logical  unit  corresponds 
in  hardware  to  the  Tektronix  screen  if  the  user  first  enters  on  the 
Tektronix,  ’CONNECT  INPUT,  OUTPUT'.  PRNPLT  is  set  FALSE  a no  CAI.PLT  is 
set  TRUE  since  plots  arc  to  be  produced  on  the  Tektronix.  TTLCP1.  is 
set  TRUE  since  a title  on  the  plots  is  to  be  made  by  using 
"SET  TITLE  = T,  X,  V,  THC.  and  LYG  will  have  their  values 

stored  every  communication  interval  through  the  PREPAR  statement.  By 
invoking  the  START  command  the  test  program  is  executed  with  the 
values  for  THDOT,  R and  THIC  as  set  in  the  INITIAL  section  of  the 
program.  The  simulation  stopping  time,  TSTP,  is  reset  to  15  sec  in  a 
run  time  SET  command.  Figures  22  and  23  show  plots  of  X and  Y and  THG 


and  LYG  ns  generated  by  the  PLOT  command.  THDOT  is  reset  to  -1.745 
rad/sec  and  t lie  program  is  run  again.  Figures  24  and  25  show  the 
resulting  plots  of  X and  Y and  THG  and  LYG.  These  tests  check  the 
operation  of  the  CPTR  macro  for  positive  and  negative  step  inputs  in 
THDOT.  The  two  D1SP1.Y  commands  result  in  data  displayed  on  the  screen 
as  well  as  recorded  on  the  'PRINT'  file.  The  PRINT  statement  causes 
recording  of  data  only  on  the  'PRINT'  file.  When  the  Tektronix  plotting 
session  is  complete  the  'PRINT'  file  is  batched  to  a printer  (not  shown 
here) . 

VIII.  CRA 

The  macro  models  a resolver  used  in  the  continuous  rotation  of 
axes  mode.  The  standard  form  is 

CRA  (X,  Y,  Til , I.Y  = V,  V,  THDOT,  THIC) 

where 

TH  = / THDOT  dt  + THIC  (as  would  be  output  from  a resolver) 

o 

X = U*COS  (Til)  - V*S 1 N(TH) 

Y = 11  * S IN  (TH)  + V*COS(TH) 

LY  = .TRllK.,  -n  TH-  n , (4n- 1 ) n ■ I'll  (4nfl)n 
( I -4n ) ii  Til  ( - I -4n)  n 
= .FALSK.,  (4n-l)n  TH  ( 4 n - 1 > u , 

( 1 — 4 n ) n I II  ( 1 -4n ) n 

n = a positive  integer. 

All  angles  and  angular  rates  are  in  radians  and  radians  per  second. 

Figure  2b  is  a flow  diagram  of  CRA  and  Figure  27  is  a listing  of 
CRA.  The  only  difference  between  this  macro  and  CPTR  is  that  RA,  the 
analog  macro  (Section  XVII)  for  flic  rotation  of  axes  is  used  instead  of 
the  ACSL  RTP  macro.  Consequently,  the  reader  is  referred  to  Section  VII 
for  a discussion  of  the  CRA  algorithm. 

Figure  28  is  the  CRA  test  program  listing.  Since  the  CPTR  test 
program  is  identical  to  the  CRA  test  program,  except  for  substituting 
U and  V for  R as  inputs  to  the  macro,  the  user  is  referred  to  Section  VII. 
Figure  29  contains  the  run  time  commands  entered  on  the  Tektronix  4014 
to  obtain  plots  and  deferred  printer  output.  These  tests  are 
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identical  to  those  described  in  Section  VII.  Figures  30  and  31  show 
plots  of  X and  Y and  THG  and  LYG  for  THDOT  input  as  a step  with  a 
magnitude  of  1.745  rad/sec.  Figures  32  and  33  contain  plots  of  X and  Y 
and  THG  and  LYG  for  a -1.745  rad/sec  step  in  THDOT.  THIC  was  0.1745  rad 
and  U and  V were  1.0  throughout  testing. 


IX.  CRTP 

The  macro  models  a continuous  rectangular  to  polar  resolver.  The 
standard  form  is 

CRTP (R , TH,  THDOT,  LY=X,  Y) 


where 


2 2 

R = X + Y 

TH  = AT AN 2 (Y,X) 

THDOT  = -r-  TH 


LY  = .TRUE.,  -h-TH-ti,  (4n-l)rr«rTH'  (4n+l)n 

(l-4n)n  TH' (-l-4n)n 
= .FALSE.,  (4n-3)n<TH- (4n-l)n, 

( 3-4n) n >TH  • ( 1 -4n) n 

n = a positive  integer. 

All  angles  and  angular  rates  are  in  radians  and  radians  per  second. 

R and  TH  are  calculated  on  both  major  and  minor  integrat ion  steps. 

THDOT  is  only  calculated  every  major  integration  step  since  some  of  the 
integration  algorithms  available  in  ACS1.  can  integrate  backward  in  time 
on  minor  integration  steps.  The  method  for  finding  THDOT  is  that  used 
in  the  ACSL  macro  DERI VT  (...).  Indeed  CRTP  calls  the  ACSI.LIR  FORTRAN 
subroutine  Z7.XDOT  to  perform  the  actual  calculation  just  as  DF.RTVT  does 
Furthermore  tin1  difference  in  I'll  between  major  integration  steps  cannot 
equal  or  exceed  3.1416  rail  with  the  algorithm  employed  in  CRTP. 


Figure  34  is  a flow  diagram  of  CRTP  and  Figure  35  is  a listing  of 
CRTP.  Referencing  the  flow  diagram,  the  ACSL  macro  RTP  converts  the 
rectangular  coordinates  X and  Y to  the  polar  coordinates  R and  THS.  If 
THS  was  changing  as  a continuous  and  monotonical ly  increasing  or 
decreasing  function,  then  it  would  look  something  like  the  curves  shown 


in  Figure  36.  Because  TH  and  LY  are  the  outputs  of  an  analog  resolver, 
they  would  be  as  shown  in  Figure  36  for  +THS.  It  is  these  variables 
that  the  remainder  of  the  algorithm  must  provide.  If  the  function 
ZZFST  is  near  zero  THDOT  is  not  calculated,  but  TH  is.  It  can  be  seen 
from  Figure  36  that  the  difference  between  TH S and  Till.,  the  value  of 
THS  at  the  end  of  the  last  major  integration  step,  is  larger  than  n rad 
only  when  there  is  a discontinuity  in  THS  (if  it  is  assumed  that  the 
rotation  of  the  point  (X,Y)  will  not  change  more  than  n rad  per  major 
integration  step).  It  X and  Y are  continuous  functions  this  assumption 
can  be  made  good  by  appropriate  choice  of  the  integration  step  size. 

The  discontinuity  in  THS  coincides  in  time  with  the  switching  of  [,Y  as 
shown  in  Figure  36.  If  the  correct  choice  of  LY  is  made  initially,  the 
value  of  LY  can  always  he  determined  by  switching  its  value  whenever 
a discontinuity  in  THS  occurs.  When  ZZFST  is  near  zero  the  difference 
between  THS  and  Till,  is  checked  (Figure  34)  to  see  it  it  is  i n 
(actually  3.1416)  rad.  II  so,  LY  is  switched  from  OLDLY , its  value  at 
the  end  of  the  last  major  integration  step.  The  macro  then  jumps  to 
the  bottom  of  the  code  where  I'll  is  set  t ii  • THS  depending  on  whether 
LY  is  TRUE  or  FALSE.  Note  that  Till,  and  OLDLY  are  not  changed  to  THS 
and  I.Y  as  would  normally  he  done  since  ZZXDOT  needs  THI.  to  be  the  value 
of  THS  at  the  end  of  the  last  major  integration  step.  Also,  OLDLY 
would  get  changed  more  than  once  per  major  integration  step  which  could 
cause  it  to  have  the  wrong,  value.  When  ZZFST  is  near  one  in  value  and 
ZZRNFL  is  TRUE,  LY  is  set  to  the  value  in  OLDLY  and  Till,  the  input  angle 
to  ZZXDOT,  is  set  to  the  value  in  TUI,.  The  macro  then  jumps  to  the 
ZZXDOT  call  to  reinitialize  ZZXDOT.  OLDLY,  THI.  and  TH  are  calculated 
as  shown  in  Figure  34.  When  ZZFST  is  near  one,  ZZRNFL  is  FALSE,  and 
Z7.ICFL  TRUE,  then  I.Y  is  set  TRUE  and  THI  is  set  equal  to  THS  and  the 
macro  jumps  to  ZZXDOT.  If  ZZICFI,  is  FAI.SF  then  THS  is  checked  for  a 
discontinuity.  If  there  is  no  discontinuity  THT  is  set  equal  to  THS 
and  the  macro  Jumps  to  the  ZZXDOT  call  to  calculate  THDOT;  if  there 
is  a discontinuity,  LY  is  switched  from  the  value  in  OLDLY  and  THI  is 
calculated  as  shown  in  Figure  36  in  order  to  preserve  the  proper  sign 
on  THDOT.  The  call  to  ZZXDOT  is  then  made  to  calculate  THDOT.  THDOTL 
and  TI.  are  not  used.  THDTIC  is  the  initial  condition  and  is  set  to 
0.0  through  a CONSTANT  statement.  OLDLY  is  reset  to  the  value  in  LY. 

THL  is  set  equal  to  THS  explicitly  even  though  ZZXDOT  sets  THI,  equal 
to  THI.  ZZXDOT  is  reliable  except  when 


THI  = '1IIS  + S 1 ON  (h.  JH  tJ,  mi.)  . 

Then  THL  will  have  the  wrong  value  on  the  next  ealeulalion  of  THDOT, 
hence  THL  is  explicitly  set  to  THS.  Til  is  set  to  'THS  depending  on 
the  value  of  LY. 

Figure  37  is  a listing  of  the  AOSI,  test  program  for  CRT!’ . In  the 
DERIVATIVE  section  the  X and  Y inputs  to  CRTI’  are  calculated  from 
Cl’TR.  Since  CPTR  has  inputs  of  R (set  to  1.0  in  the  INITIAL  section) 


and  THD,  the  angular  rate,  then  CRTP  should  output  values  in  R and 
THDOT  that  are  similar  to  1.0  and  THD.  Similarly,  THS,  the  angular 
output  of  CPTR,  should  match  TH.  THD  is  a step  function  of  magnitude 
K (radians)  if  LSTEP  is  TRUE  and  is  a sine  wave  of  magnitude  K (radians) 
and  frequency  W (rad/sec)  if  LSTEP  is  FALSE.  TH1C  as  defined  in  the 
INITIAL  section  is  the  initial  condition  (radians)  for  the  CPTR  inte- 
gration of  THD.  LYG  has  a value  of  1.0  if  LY  is  TRUE  and  a value  of 
0.0  if  LY  is  FALSE.  THDG , THSC. , THG  and  THDOTG  are  THD,  THS,  TH,  and 
THDOT  expressed  in  degrees  or  degrees  per  second.  The  run  time  commands 
entered  on  the  Tektronix  for  hard  copy  plots  and  printer  output  of  the 
test  results  are  shown  in  Figure  38.  LSTEP  is  set  FALSE  to  input  THD 
as  a sine  wave.  THIC  is  set  to  -14.31,  the  appropriate  initial  condition 
for  a cosine  wave  of  magnitude  K/W  where  K is  set  to  1.431  rad  and  W 
remains  0.1  rad/sec.  TSTP,  the  simulation  stopping  time,  is  set  to 
70  sec  to  ensure  a full  cycle  of  the  THD  sine  wave.  The  input  sine 
wave  tests  the  CRTP  logic  for  angle  changes  due  to  X and  Y that  are  not 
monotonical ly  increasing  or  decreasing.  The  program  is  executed 
through  the  START  command.  The  PLOT  command  is  used  to  generate  plots 
of  THG,  THDOTG,  LYG,  THDG  and  THSG  as  shown  in  Figures  39,  40  and  41. 
Figures  39  and  41  are  identical  in  THSG  and  THG  and  are  approximately 
the  same  in  THDG  and  THDOTG.  The  run  time  command 

PRINT  'ALL' 

records  data  specified  in  the  PREPAR  statement  onto  the  'PRINT'  file 
which  can  be  batched  to  a line  printer  later. 

K is  reset  to  1.745  and  THIC  is  reset  to  -17.45  and  the  program 
is  run  again.  Figures  42,  43  and  44  show  plot  results  for  THG,  THDOTG, 
LYG,  THDG  and  THSG.  The  basic  difference  between  this  test  and  the 
last  is  the  starting  quadrant  of  THG.  These  tests  were  made  to  ensure 
CPTR  and  CRTP  handle  different  starting  quadrants  correctly. 

LSTEP  is  set  TRUE,  THIC  is  set  to  0.0  and  TSTP  is  set  to  20.0  sec 

to  run  a step  function  in  THD  for  a test  with  angle  changes  due  to  X 

and  Y increasing  monotonical 1 y . Plots  resulting  from  these  tests  for 
THG,  THDOTG.  X,  Y,  THSG  and  LYG  are  shown  in  Figures  45,  46  and  47. 

K is  set  to  -1.745  to  perform  tests  with  monoton ical ly  decreasing  angle 
changes  due  to  X and  Y.  Figures  48,  49  and  50  are  plots  of  THG,  THDOTG, 
X,  Y,  THSG  and  LY  for  these  tests. 

K is  set  to  1.5708  and  -1.5708  to  test  for  any  changes  in  TH  due 

to  it  being  very  near  multiples  of  n(3.1416)  rad,  that  is,  the  points 

where  LY  is  being  switched.  No  occurrences  of  this  sort  happened  as 
can  be  seen  in  Figures  51  and  52,  which  are  plots  of  THG  and  THDOTG  for 
the  two  different  values  of  K. 
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The  operation  of  CRTP  through  a REIN1T  run  time  command  was  tested 
by  running  the  program  with  TSTP  set  at  2.2  sec,  invoking  a REINIT, 
resetting  TSTP  to  6.1  and  again  executing  the  program.  Both  runs  were 
later  printed  and  compared  tor  the  correct  transition  through  the  time 
REINIT  was  invoked.  This  test  showed  that  values  immediately  before 
and  after  the  REINIT  command  were  the  same  and  therefore  the  test  was 
successful . 


X.  DIFF 

The  macro  models  a logic  differentiator.  The  standard  form  is 

LY  = DIFF (LX). 

LY  is  a logical  variable.  l.X  is  any  logical  expression.  I.Y  is  TRl'E  for 
one  calculation  interval  (clock  period  on  the  analog)  following  the  time 
that  LX  changes  f rom  FALSE  to  TRUE.  I.Y  does  not  become  TRUE  again  while 

LX  remains  TRUE.  It  is  to  be  noted  that  this  is  a leading  edge  differ- 

entiator only. 

Figure  53  is  a t low  diagram  ot  DIFF  and  Figure  34  is  a I ist ing  of 
DIFF.  Figure  5 3 is  believed  to  he  self-explanatory  except  when  7.7. 1 CFL 
is  TRUE  and  7.ZRNFI.  is  FALSE.  Obviously  LY  should  he  FALSE  for  this 
case,  but  LX I.  had  to  he  made  FALSE  also  to  ensure  that  LY  went  TRUE  for 
the  case  of  LX  being  TRUE  at  zero  t imo. 

Figure  55  contains  a I ist  ing  ot  the  DIFF  At'SL  test  program  and  run 
time  commands.  LX  is  made  I'RUF  or  FALSE  as  a fund  ion  ot  time  bv 
alternately  adding  and  subtracting  SITE  macros  from  one  another  and 
checking  to  see  if  the  results  are  approx  imat e 1 v either  one  or  zero. 
Instead  of  inserting  I.X  as  the  input  argument  in  DIFF,  its  equivalent 

expression  was  inserted  to  help  check  lor  proper  macro  operation.  NDBUG 

is  set  to  8 to  inspect  the  oper.il  ion  o|  DIFF  through  intermediate 
integration  steps.  TSTP,  the  simulat  ion  stopping  t imo,  is  set  to 
6.0  sec.  The  simulation  is  run  bv  invoking,  the  START  command.  REINIT 
is  commanded,  TSTP  is  reset  to  10. ti  so  and  the  program  is  run  again. 
Both  runs  are  printed  e it  ..coining  to  the  l’REPAR  statement  as  shown 
i:i  Figure  56. 

XI.  DLYFF 


The  macro  models  a delay  I I ip  Mop.  The  standard  form  is 
LY  = DI.YFF(LX). 

LY  is  a logical  variable  and  LX  is  any  logical  expression.  LY  follows 
the  value  of  I.X  with  a delav  ol  one  calculat  ion  interval  (clock  period 
on  the  analog). 


Figure  57  is  a flow  diagram  for  DLYFF  and  Figure  58  is  the  macro 
listing  for  DLYFF.  Figure  59  is  the  test  program  listing  with  run  time 
commands.  In  the  test  program  DERIVATIVE  section  LX  is  determined  by 
whether  a function  composed  of  STEP  macros  is  near  zero  or  one.  In  the 
run  time  command  section  NDBUG  is  set  to  8 and  the  program  is  run  twice 
with  a REINIT  statement  between  runs.  The  resulting  output,  in  accord- 
ance with  the  PREPAR  statement,  is  shown  in  Figure  60.  It  should  be 
noted  that  the  one  calculation  interval  lag  seen  in  l.Y  appears  to  be 
two  intervals.  This  is  due  to  LY  being  printed  TRUE  at  the  end  of  3 sec, 
even  though  it  was  TRUE  between  2.0  and  3.0  sec  as  the  NDBUG  printout 
of  the  first  8 integration  steps  would  show. 


XII.  LMMINT 

The  macro  models  a limited  mode  controlled  integrator.  The 
standard  form  is 

y = LMMINT  (YD,  IC,  LI,  LO,  LI.,  UI.) 

where  YD  is  the  derivative  expression  and 

IC  = Y(o)  - same  restriction  on  IC  as  the  INTEG  statement 

LI  and  LO  are  logical  variables  or  expressions  of  arbitrary  complexity 
denoting  the  mode.  The  truth  table  is 


LI 

LO 

MODE 

T 

F 

RESET 

F 

F 

OPERATE 

where  T = TRUE 

T 

T 

OPERATE 

and  F = FALSE 

F 

T 

HOLD 

LL  is  the  lower  limit  on  Y and  UL  is  the  upper  limit  on  Y. 

Figure  61  is  a flow  diagram  for  LMMINT  and  Figure  62  is  the  macro 
listing  for  LMMINT.  The  following  comments  parallel  the  flow  diagram. 
YDOT  is  set  equal  to  YD  so  if  YD  is  an  expression  it  will  be  calculated 
only  once.  If  ZZFST  is  not  less  than  0.5  and  if  ZZRNFI,  and  ZZICFL  are 
FALSE,  then  Y is  checked  for  a value  that  lies  outside  (l.L,  UL) . When 
Y does  lie  outside  (LL,  UL)  then  MODE,  the  logical  variable  used  to 
zero  the  derivative  of  STATE,  is  calculated.  If  the  derivative  YDOT  is 
of  such  a sign  as  to  drive  Y further  into  the  limit  then  MODE  is  set 
TRUE.  With  MODE  TRUE;  the  derivative  of  STATE, as  calculated  by  the 
ACSI.  RSW  function,  is  zero  to  inhibit  driving  Y further  into  the  limit. 
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When  YDOT  is  of  such  a sign  as  to  cause  Y to  come  off  a limit,  MODE  will 
he  FALSE,  thus  causing  the  STATE  derivative  to  be  set  to  YDOT.  M0DE1 , 
the  reset  mode  logical  variable,  is  calculated.  If  M0DE1,  MODE  or  the 
logical  expression  ( (I.T)  .AND.  .NOT.  (1.1))  is  TRUE  then  the  derivative  of 
STATE  is  set  to  0.0  by  RSW.  11  MODE1  is  TRUE,  l CPF,  the  variable  used 
to  bias  Y back  to  its  initial  value  ICP,  is  set  to  STATE.  Y is  calculated 
using  STATE,  ICP,  and  1CPP.  When  ZZFST(YDP)  is  less  than  0.5  the  macro 
calculates  STATE  using  the  previous  RSW  and  Y is  Evaluated.  It  is  to  be 
noted  that  when  Z.ZFST  is  less  than  0.5  the'  macro  uses  whatever  value 
for  MODE  that  was  calculated  in  the  first  intermediate  integration  of 
this  major  integration  step  to  determine  the  derivative  of  STATE,  no 
matter  what  sign  changes  have  occurred  in  YD  during  the  interim.  If 
ZZRNFI.  is  TRUE  ICP  is  set  to  Y,  1CPP  is  set  to  STATE,  and  FLAG  is  set 
to  1.  If  ZZ1CF1.  is  TRUE  then  MODE  is  set  FALSE,  and  if  FLAG  has  a 
value  of  0 (which  it  will  for  the  first  pass  through  the  INITIAL 
section)  then  ICP  is  calculated.  Using  ICP  instead  of  1C  forces  the 
initial  value  ot  Y to  lie  in  the  open  interval  (LL,  UL)  (A  limited 
analog  integrator  will  not  output  values  above  or  below  its  limits 
regardless  of  the  value  to  which  the  l.C.  pot  is  set.).  The  macro 
then  branches  to  calculate  Y and  STATE. 

Figure  tit  is  tin1  macro  test  program  list  ing.  The  output  derivative 
to  LMM1NT  in  t hi'  DERIVATIVE  sect  ion  is  a sine  wave  with  magnitude  A and 
frequency  I'W0P1*WCPS  rad/sec  as  defined  in  the  INITIAL  section.  Since 
Y will  be  a cosine  wave  of  magnitude  A/ (TWOP1 *WCPS) , the  1C  is  set 
according.lv  in  the  INI  II  \l  section.  1.1.  and  UL  are  defined  in  the 
INITIAL  sect  ion  as  -2.  ind  F2.  1, 1C  and  I.TG  cause  l.l  and  LT  to  be  TRUE 
or  FALSE  depending  on  il  t hoy  arc  1.0  or  0.0.  LIG,  I.TG  and  YD  are 
used  to  represent  II,  I.T  and  the  I .MM  l NT  derivat  i ve  argument  for  plotting 
purposes.  Figure  li-<  is  a listing  ot  the  run  t into  commands  entered  on 
t ho  Tektronix  401 A to  obtain  hard  copv  plots  ol  the  results,  as  well 
as  generating  a 'PRINT'  I i I < to  be  later  output  on  a line  printer. 

A series  of  REINIT  commands  weio  used  to  check  the  REINIT  actions 
under  the  limit  , hold  and  reset  modes.  In  order  to  obtain  plots  on 
the  Tektronix,  CAI.I'I.I  is  set  TRUE  and  PRNIM.T  is  set  FALSE.  To  get 
plots  with  titles,  TTLCPI  is  set  TRUE  and 

SET  TITLE  = ’ . . . ' 

is  invoked.  I'S'TP,  the  sitniilat  ion  slopping  t into  was  sot  to  0.44  sec; 
the  simulation  is  run  via  the  START  command  and  stopped  via  the  TF.RMT 
macro  in  the  tost  program  DYNAMIC  sect  ion.  The 

SAVE  'BASE' 

command  saves  the  init  ini  run  rend i I ions  (otherwise  they  would  be  lost 
after  a REINIT  command)  so  that  thov  can  be  restored  later  for  other 
tests.  REINIT  is  invoked  and  i'STP  is  reset  to  1.6  sec.  Instead  of 
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plotting  and  printing  these  runs  separately  NRWITG  is  set  TRUE  so  that 
the  file  containing  the  values  of  all  the  variables  on  the  PREF  >.  list 
is  not  rewound  between  runs,  thereby  accumulating  all  of  the  data  from 
the  runs  as  a unit  to  be  plotted  and  printed  together.  The  ACTION 
command  causes  NDBUG  to  be  set  to  8 when  time  is  S 2.191  sec.  After 
TSTP  is  set  to  3.6  sec  the  ACTION  is  cleared.  The  simulation  variables 
are  printed  for  the  next  eight  integration  steps  to  see  exactly  when 
LIG  and  LTG  switch  their  values  to  ensure  that  the  LMMINT  macro  is 
changing  modes  at  the  right  point  in  time.  The  simulation  is  run  by 
the  START  command  followed  by  a REINIT  command  for  TSTP  values  of 
0.44,  1.6  2.9  and  3.6  sec  to  catch  the  macro  in  its  various  modes  with 
REINIT  invoked.  The  last  run  for  this  series  of  tests  starts  at  3.6  sec 
and  ends  at  5.0  sec.  NRWITG  is  set  FALSE  and  Y,  YD,  LIG  and  LTG  are 
plotted  with  the  PLOT  command  as  shown  in  Figures  65,  66  and  67.  The 
PRINT  command  causes  the  variables  in  the  PREPAR  statement  to  be 
written  to  the  'PRINT'  file  which  is  later  batched  ti  a line  printer 
for  listing.  The  proper  operation  of  LMMINT  under  REINIT  conditions  is 
checked  by  simplv  saying  START  without  an  accompanying  REINIT.  TSTP  is 
set  to  5.1  sec  and  the  data  are  printed.  (No  data  are  shown  in  this 
document.)  The  initial  conditions  were  correctly  established  for  Y at 
the  simulation  starting  time  of  3.6  sec  (the  last  Lime  a REINIT  command 
was  invoked . ) . The 

RF.STOR  'BASE' 

statement  restores  in  the  simulation  the  original  initial  conditions. 

A is  set  to  6.2832  so  that  the  output  Y does  not  limit  for  this  run. 

TSTP  is  set  to  5.0  and  START  commences  running  the  simulation.  PLOT 
commands  that  a graph  of  Y and  YD  he  made.  The  result  is  shown  in 
Figure  68.  LIG  and  LTG  are  as  before  and  are  therefore  not  shown. 

XIII.  LPULSE 

A train  of  logical  pulses  can  be  generated  using  the  LPULSE 
macro.  The  independent  variable,  default  T,  is  used  to  drive  it.  It 
is  to  be  noted  that  in  an  all  digital  CSSI.  implementation  and  in  digital 
regions  in  a hybrid  implementation  the  integration  step  size  may  affect 
the  answers  in  that  a too  large  step  could  cause  the  pulse  to  stay  on 
indefinitely.  The  output  will  always  be  high  (equal  to  TRUE)  at  the 
beginning  of  the  first  calculation  interval  (clock  period  on  the  analog) 
that  follows  the  exact  turn-on  time.  The  standard  form  is 

LY  = LPULSE  (tz,  p,  w) 

The  result,  LY,  is  a logical  pulse  train  (TRUE  or  FALSE)  starting  at 
the  first  calculation  interval  that  equals  or  exceeds  tz.  Period  is 
p and  width  is  w.  Figure  69  is  the  macro  listing  of  LPULSE.  Figure  16 


is  the  test  program  listing  tor  I. PULSE . I.Y  2 is  the  output  from  LPULSE 
us  invoked  in  the  DERIVATIVE  section.  The  variables  tz,  p and  w are 
set  by  a CONSTANT  statement  in  the  INITIAL  sect  ion  to  2.0,  4.0  and 
1.0  see,  respectively.  The  results  of  the  test  are  shown  in  Figure  17. 
Section  VI  describes  the  run  t ime  commands  for  the  test . 

XIV.  LSTEP 

The  LSTEP  macro  produces  a change  f i om  FALSE  to  TRUE  in  the 
output  at  a specified  value  ol  the  independent  variable  T.  The  standard 
form  is 

I.Y  = LSTEP ( t z) 

where 


I.Y 

= FALSE  T t 

I.Y 

= TRUE  r t 

it* 

log 

ic  step  start s 

.a  tlu  1 i rs 

t i a 1 ni  1 a t i 

.’il  into 

rva  1 

(clock  period 

on  the  analog)  that  equals  or  exceeds  t /. . 

Figure  70  is  the  macio  I ist  ing  ol  LPULSE.  Figure  1 1>  is  the  test 
program  listing  tor  I.SIT'.I’.  As  shown  in  the  DERIVATIVE  section  LYf  is 
the  output  c)l  a logical  expression  involving,  on  I \ the  LSTEP  macro.  Two 
runs  are  made  with  the  test  program,  the  second  run  start  ing  with 
a REINIT  command.  (Sect  ion  VI  describes  the  run  time  commands.)  The 
result  is  shown  in  Figure  I/. 

XV.  MODE 

Reference  to  the  ana  I op,/ I og  i c modes  is  bv  means  ol  this  pseudo 
simulation  operator.  The  standard  term  is 

I.Y  = MODE  (a  rg) 

It  has  no  meaning,  in  an  all  digit. il  OSSI.  implemenlat  ion  and  is  always 
FALSE.  For  analog  imp  lcment.it  ion,  the  operator  is  TRUE  when  the 
named  mode  is  in  oper.it  ion. 

arg  - CI.R , STP , RUN  ,S  ,M , E .MS  , SEC  , K'l  , ST,  PS  , I ('  ,M  or  OP 

Figure  71  is  a I ist  ing  ol  the  macro  MODE.  Figure  It)  is  the  test 
program  listing.  I.Y  1 is  the  output  variable  from  MODE,  as  invoked 
in  the  DERIVATIVE  sect  ion.  X is  the  dummy  input  variable  to  MODE. 
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Two  runs  are  made  with  the  test  program,  the  second  starting  with  a 
REINIT  command  (Section  VI  describes  the  run  time  commands.)  The 
result  is  shown  in  Figure  17. 


XVI.  MONO 

The  macro  models  a logic  monostable.  The  standard  is 

LY  = MONO(R.l.X) 

where  l.Y  is  a logical  variable;  IX  is  i variable,  logical  expression, 
or  constant;  and  R is  a real  constant,  expression  or  variable  defining 
the  period  of  t he  monostable.  l.Y  is  TRt'K  lor  a period  of  R seconds 
following  the  time  that  I.X  becomes  TRUE.  After  this  period  l.Y  becomes 
FALSE  unless  I.X  is  still  TRUE.  If  I.X  goes  FALSE  during  t ho  period, 

LY  still  remains  TRUE  for  the  rest  ol  the  period.  Figure  72  is  a flow 
diagram  for  MONO  and  Figure  71  is  a list  ing  of  t hi*  macro.  Following 
the  flow  diagram  to  where  ZZ.  ICEI.  is  TRIM.,  l.Y  is  set  equal  to  LX.  If 
LY  is  TRUE  GLUT,  the  value  ot  1 when  LY  lirst  went  lKl'K  in  this  period, 
is  set  to  T.  (NOTE:  It  is  implicit lv  assumed  in  this  macro  that  the 
independent  variable  has  its  default  name,  'T'.)  Iletoic  leaving  the 
macro  the  ACSL  Executive  branches  t o set  01.010,  the  value  for  the 
previous  LY,  to  l.Y.  When  ZZICEI.  is  FALSE  EXE  is  set  equal  to  LX  so  if 
LX  is  an  expression  it  will  have  to  be  evaluated  only  once  per  macro 
call.  If  OLDLY  is  FALSE  I.XI*  is  checked.  If  EXE  is  FALSE,  LY  is  set 
to  FALSE,  OLDLY  is  sot  equal  to  l.Y  and  the  macro  is  terminated.  If  LXE 
is  TRUE,  OI.DT  is  set  to  T,  l.Y  is  set  to  TRUE,  OI.DI.Y  is  sot  to  l.Y  and  the 
macro  is  terminated.  If  OI.DI.Y  is  TRUE,  I is  subtracted  from  OI.DT  to  see 
if  the  period,  R,  has  been  exceeded.  II  the  period  has  been  exceeded, 

LXE  is  checked  and  branching,  occurs  as  beloro.  It  the,'  period  has  not 
been  exceeded,  LY  is  set  to  I'RUE,  OI.DI.Y  is  set  to  l Y and  the  macro  is 
terminated.  If  ZZEST(OI.DI.Y)  is  loss  than  O.S  or  ZZ.RNFI.  is  TRUE, the 
macro  code  is  skipped. 

Figure  74  is  a list  ing  of  the  test  program  and  associated  run  time 
drive  cards  tor  MONO.  I.X,  the  logical  input  to  MONO , is  calculated  in 
the  DERIVATIVE  section  where  an  expression  of  ACSL  STEE  operators  being 
greater  than  O.1)  or  not  determines  whether  I.X  is  TRUE  or  FALSE.  In  the 
INITIAL  sect  ion,  R is  set  to  2.0  sec  through  the  use.'  of  the  CONSTANT 
statement.  The  MONO  is  invoked  in  the  DERIVATIVE  sect  ion  with  LY  as 
its  output. 

In  the  run  time  commands  I'STE,  the  simulation  stopping  time,  is 
set  to  7.0  sec.  The  simulation  is  started  with  the  START  command  and 
is  terminated  by  the  TKRMT  macro  located  in  the  program  DYNAMIC  section. 

The  ER1NT  command  lists  the  variables  as  dictated  bv  the  PREEAR  statement. 
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The  REINIT  command  is  invoked,  TSTP  is  reset  to  12.0  sec,  the  simulation 
is  executed  and  the  variables  in  the  PREPAR  statement  are  printed  again. 

Figure  75  is  the  data  output  from  the  line  printer.  The  monostabie 
period  and  logic  output  are  correct  and  the  macro  performed  as  it  should 
under  the  REINIT  condition. 

1 

XVII.  RA 

The  macro  performs  a rotation  of  axes.  The  standard  form  is 
RA(X,Y»U,  V.TII) 

where 

if 

' x = u*cos  mo  - v * s i n ( nn 

Y U*S  IN  (Til)  + V*('OS(TH) 
ill  = the  angle  in  radians 

Figure  7b  is  a listing  ot  the  macro.  This  macro  was  checked  out 
simultaneously  with  the  ERA  macro  as  previously  described  in  Section  VIII. 

XVIII.  SHIFT 

The  macro  models  a '(-hit  shill  reg.ister.  The  standard  form  is 

| 

Sll  I FT(  1.Y0,  IV  1 .I.Y.’.I  Y ! NIC,  I.XSII , I.XS  I ) 

where  1.Y0 , 1.Y  1 , 1.Y2 , l,Y  I are  the  logical  output  variables  corresponding  to 
the  lour  flip-flops  in  the  register  and  represent  I , 2 , 4 , 8, respec t i ve 1 y . 

N I C is  an  integer  (0  to  15)  corresponding  to  a hit  pattern  which 
represents  the  register  initial  state.  I.XSII  and  I.XSI  are  logical 
variables  or  logical  expressions.  Tin  register  shifts  its  hit  pattern 
to  the  left  (I.YO  I.XSI.  I.YI  I.YO,  I.Y2-I.YI,  I.Y  1 I.Y2)  every  major  inte- 
gration step  during  which  I.XSII  is  TRUE,  otherwise  the  present  state  is 
held. 

Figure  77  contains  a I low  diagram  of  SHIFT.  Figure  78  is  an  ACSI. 
listing  of  SHIFT.  In  the  I low  diagram  it  X.Z.FST(K)  is  less  than  0.5  or 
ZZ.RNF1.  is  TRUE  then  the  macro  code  is  bypassed.  Hence  this  macro  is 
executed  only  once  per  major  integral  ion  step  and  when  the  REINIT  run 
time  command  is  invoked  no  values  are  changed.  II  ZZ.ICFl.  is  TRUE, 

I.YO  through  I.Y3  are  initialized  by  AND'ing  NIC  with  1,2, 4, 8,  respect- 
ively. It  Z.Z.IC.FI,  is  FAI.SE  and  I.XSII  is  FAI.SF  the  rest  of  the  macro  code 
is  skipped.  IF  I.XSII  is  TRUE  the  register  is  shifted  to  the  left  bv 
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one  bit.  LLYO  through  1.LY3,  the  old  values  for  the  register  before  a 
shift  occurred,  are  set  to  LYO  through  1.Y3,  respectively. 

Figure  79  is  a listing  of  the  ACSI.  test  program  for  SHIFT.  In  the 
DERIVATIVE  section  I.XSl  and  l.XSH  are  made  functions  of  time  using  the 
STEP  macro.  NIC  is  set  equal  to  15  in  the  INITIAL  section  via  a 
CONSTANT  statement.  SHIFT  is  coded  in  the  DERIVATIVE  section  with  the 
same  variable  names  as  used  above. 

In  the  run  time  commands  the  variables  to  be  printed  are  listed 
in  a PREPAR  statement.  TSTP,  the  simulation  stopping  time,  is  set  to 
5 sec.  The  simulation  is  started  and  the  resulting  data  printed  in 
accordance  with  the  PREPAR  statement.  The  REIN1T  command  is  invoked, 
TSTP  is  set  to  11  sec  and  the  simulation  is  run  again  with  a printout 
of  data  as  before.  Figure  80  gives  tin-  results  as  output  on  a line 
printer. 


XIX.  SRTEFF 

The  macro  models  a flip-flop  in  its  set,  reset,  trigger  and 
enable  modes.  The  standard  form  is 

LY  = SRTEFFf  LI  C , LS  , I.R , LE) 

L1C  is  the  logic  expression  for  the  init ial  state  ol  the  flip-flop. 

The  flip-flop  changes  state  only  when  the  logic  expression  LE  is  TRUE. 
The  logic  variable  LY  is  controlled  by  the  logic  expressions  LS  and 
I.R  as  shown  in  the  following  truth  table. 


LS  I.R  I.Y 

F F T/F  1.Y  is  unchanged 

T F T 

F T F 

T I T/F  1A  changes  state 


Figure  81  contains  a I low  diagram  ot  SRTF.FF.  Figure  82  is  an  ACSI. 
listing  of  SRTF.FF.  Referencing  the  flow  diagram,  if  ZZFST(OLDLY)  is 
less  than  0.5  or  77.KNF1.  is  TRUE  the  macro  code  is  bypassed.  Thus  this 
macro  is  executed  only  once  per  major  integration  step.  If  7.7. 1 CFL  is 
TRUE,LY  is  set  equal  to  Lit'  and  control  is  transferred  to  the  last 
executable  statement  in  the  macro.  II  1.1  is  FALSE  the  rest  of  tin*  code 
is  skipped.  The  rest  of  the  branching  straightforwardly  implements  the 
truth  table  for  LY . However  there  is  the  subtlety  of  the  logic  sequence 
being  important.  Finally  OI.DY,  the  previous  value  of  Y,  is  set  equal 
to  I.Y  in  preparation  for  the  next  pass  through. 
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Figure  83  is  a listing  of  the  ACSL  test  program  for  SRTEFF.  In 
the  DERIVATIVE  section  LE,  LR,  LS  are  made  functions  of  time  via  the 
STEP  macro.  SRTEFF  is  then  invoked.  1. IC  is  defined  in  the  INITIAL 
section  through  the  use  of  the  CONSTANT  statement. 

In  the  run  time  commands  the  variables  to  be  printed  are  listed  in 
a PREPAR  statement.  TSTP,  the  simulation  stopping  lime,  is  set  to 
9 sec.  The  START  command  causes  the  simulation  to  begin  and  the  result- 
ing data  printed  in  accordance  with  the  PREPAR  statement.  The  REIN1T 
command  is  then  used  with  TSTP  being  reset  to  14  sec.  The  simulation 
is  run  again  with  a printout  of  data  as  before.  Figure  84  contains  the 
data  printout  for  these  two  runs. 
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XX.  SRTFF 

The  macro  models  a flip-flop  in  its  set,  reset  and  trigger  modes 
(no  enable  control).  The  standard  form  is 

LY  = SRTFF( I.IC,  I.S  , I.R) 

l.IC  is  a logic  expression  lor  the  initial  state  ol  the  flip-flop.  The 
logical  variable  I.Y  is  controlled  by  the  logic  expressions  LS  and  LR  as 
shown  in  the  following  truth  table. 


LS 

LR 

LY 

F 

T 

F 

T/F 

T 

I.Y  is  unchanged 

1 

F 

T 

1 

F 

T 

T 

T/F 

1.)  changes  state 

The  flow  diagram  tor  SRTFF  is  the  same  as  that  for  SRTEFF, 

Figure  81,  except  that  there  is  no  cheek  foi  the  st.ite  of  the  variable 
LE.  Figure  83  is  an  ACSL  1 ist  ing  ol  SRTFF.  The  test  program  for 
SRTEFF  is  the  same  one  used  lor  SRTFF  (Figure  St).  For  an  explanation 
.’I  the  macro  codi  anil  the  test  program  one  is  referred  to  Section  XIX. 
Figure  8b  contains  the  test  data  listing  for  SRTFF. 

XXI.  SWITCH 

The  macro  models  the  switching  lor  a single  input/single  output 
divice.  The  standard  form  is 

Y = SWITCH  (LX,  X) 
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where 


Y = X If  LX  is  TRUE 

Y = 0.0  if  LX  is  FALSE 

and  Y and  X are  taken  to  be  a real  variable  and  a real  expression, 
respectively. 

Figure  87  is  the  ACSL  code  for  SWITCH.  The  RSW  macro  is  used  with 
one  argument  being  0.0.  The  test  program  is  shown  in  Figure  16.  INPUT, 
the  switched  input,  and  LX  are  defined  in  the  DERIVATIVE  section. 

SWITCH  is  invoked  with  the  output  being  Y4 . The  run  time  commands  set 
constants,  print  and  plot  data  and  control  execution  of  the  simulation. 
The  PRF.PAR  statement  contains  the  variables  to  be  printed  as  a function 
of  time.  NORUG  is  set  to  8.  This  causes  the  printing  of  most  of  the 
program  variables  for  tin1  first  eight  integration  steps.  TSTP,  the 
simulation  stopping  time,  is  set  to  2 sec.  The  START  command  causes 
the  simulation  to  begin  and  the  result  ing  data  are  printed  in  accordance 
with  the  PRF.PAR  statement.  The  RF.INIT  command  is  then  used  with  TSTP 
being  reset  to  8 sec.  The  simulation  is  executed  again  and  data  are 
printed.  Figure  17  shows  Y4,  LX  and  INPUT  listed  as  functions  of  time. 


XXII.  TIMER 

The  macro  models  an  interval  timer  (logical  counter).  The 
standard  form  is 

TIMER  (LY , I.YCO  = NIC,  1.XC1  , IS.  I.K) 

NIC  is  an  integer  expression  spec  il’v  ing  the  number  of  pulses  that  will 
be  counted.  LS , a logic  expression,  starts  the  timer  off.  LY,  the 
logic  output  variable,  will  then  go  TRUE  and  the  count  starting  with  a 
value  equal  to  NIC  is  decremented  each  t ime  t ho  logic  expression  LXC1 
is  TRUE  on  the  first  minor  step  of  a major  integration  step.  LY  returns 
to  FALSE  the  calculat  ion  interval  after  the  one  where  the  count  reaches 
zero.  LYCO,  a logic  output  variable,  is  TRUE  for  one  major  integration 
step  (clock  period  on  the  analog)  prior  to  IY  returning  to  FALSE,  that 
is,  LYCO  is  TRUE  when  t hr-  count  is  zero.  Anv  time  that  the  logic 
expression  LR  is  TRUE,  LY  returns  to  FALSE  and  the  counter  is  reset  to 
NIC. 


Figure  88  contains  a Mow  diagram  of 
listing  of  TIMER.  Referencing  the  flow  d 
than  0.5  or  Z7.RNFL  is  TRUE  the  macro  code 
executed  only  once  per  major  integral  ion 
command  is  invoked  no  values  are  changed. 


TIMER.  Figure  89  is  an  ACSL 
iagram,  il  7.ZFST  (KP)  is  less 
is  bypassed.  This  macro  is 
step  and  when  the  REM  NIT 

If  7.7'ICFL  or  LR  is  TRUE  then 
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LY  and  STARTD,  the  local  start  down  count  logic  variables  are  set  to 
FALSE.  Also  KP,  the  local  integer  variable  that  keeps  track  of  the 
previous  value  for  the  down  count,  is  set  equal  to  NIC.  If  STARTD  and 
LS  are  FALSE,  branching  is  as  above.  If  either  STARTD  or  LS  is  TRUE, 
STARTD  and  LY  are  set  to  TRUE.  If  LXCI  is  FALSE  the  remainder  of  the 
macro  is  skipped;  if  TRUE  K,  the  present  value  for  the  down  count,  is 
decremented  by  one.  If  K is  equal  to  zero  then  LYCO  is  set  TRUE, 

STARTD  is  set  FALSE  and  KP  is  reset  to  NIC;  if  K is  not  equal  to  zero 
KP  is  set  to  K and  the  macro  is  ended.  LYCO  stays  TRUE  for  only  one 
calculation  interval  by  setting  it  to  FALSE  at  the  front  of  the  macro. 

I.S  has  to  be  TRUE  for  onlv  one  calculation  interval  for  down  counting 
since  STARTD  will  be  set  TRUE  until  K becomes  zero. 

Figure  90  is  a listing  of  the  ACSI.  test  program  for  TIMER.  In  the 
DERIVATIVE  section  l.XCl,  LS  and  LR  are  made  functions  of  time  through 
the  use  of  the  STEP  macro.  In  tin*  INITIAL  section  NIC  is  declared  an 
integer  and  set  to  7 through  the  use  of  a CONSTANT  statement.  The 
TIMER  macro  is  also  put  in  the  DERIVATIVE  section. 

In  the  run  time  commands  the  variables  to  be  printed  are  listed  in 
a PREPAR  statement.  TS'TP,  t hi'  simulation  stopping  time,  is  set  to 
18  sec.  Thi-  simulation  is  started  and  the  resulting  data  printed  in 
accordance  with  the  PREPAR  statement.  The  RE  I NIT  command  is  invoked 
and  TSTP  is  reset  to  111.9  sec.  NDRUC  is  set  to  r>  and  when  the  START 
command  is  again  invoked  a listing  of  most  ol  the  simulation  variables 
will  be  made  tor  the  tirst  r>  integration  steps  ot  this  simulation  run. 
Data  are  also  printed  according  to  the  PREPAR  statement  by  again  invoking 
the  PRINT  statement.  Figure  91  shows  the  line  printer  listing  for  these 
two  simulation  runs  (minus  the  NDRUC  portion!. 
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p . I Ki  R i T 

ATTACH. KACriL  .DC^ACF  l!  . TOrftMCM  CYS 
ATTAC*-‘»ACcl.  .DC Af SI  . j'.'TOr acsi 
ACSL  » I = l * PUT 

RFTUPM. ATCl  .maCT  II 
FTK,  I=fOMHIU'.P  = ? 

A T T AC**  • AC  CLL  IP.UCAC  SI  i.  !•*.  IO  = nC-'Cc|  SYS 
XFQ 

I.PSF  T .1.  1 P-Af  SLL  Ir>.P^*-cp  1 r Tt.PFT 

LOAO  = l.r-C 

MOGO=l  POP 

PFTI’P'y* AC^Ll  I p 

Lr  oh 

Figure  1.  NOS/BE  Interactive  terminal  commands  to 
compile  and  execute  a macro  test  program. 


P . I NPI  >T 

ATTACH. PACE  II  .DCVACFIl  .1  > = hr/,CM_SY'-> 

ATTACH.  ACFL.hC  ACSI  . TU  = '>r  Af^l  SY<- 
AC  SI  • 1 - INPUT 
PF  Tl'PU.ACSL  . MACf  II 
FTN. I=C0VPILF • *v - A 

ATTACH. ACSLl  lP.riCACFll  I •-  . | 1 ;-PC  Af"r  I Ws 
aTTACH.PL^Tt  I'J.  Tr  KTP  I * .0  1 U • I h=  ./  T H»l  ()T  ,(  y = ? 
PFOI  'F.  ST  .1  C 0M  . »*'f 
X F 0 

l rSFT  .1  IPrPL  OTL  I p .S'Mct  - //OF  A^-Tr  Khl  T 
l CST  T.L  IP-ACbl  LlH.PP*  ct  T-TSrFT 

l OAO=I.PO 
t"0C,0=l  COP 

C A T AL  Of- , l*'  0“  » k‘  A.CpO  tc<-  tc  , | p = hn  \ y x 1 1 

PF  TiiPM.PL  OTl  I p . ACSL  I I°.l  '-0" 


Figure  2.  NOS/BE  interactive  terminal  commands  to 
obtain  an  absolute  file,  t.COB,  of  a macro  test 
program  to  run  tests  for  a macro  on  tbe  Tektronix 
4014 . 
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0000000000000000000000 
HOMOO  « CM 77000* 

ACCT <PN»MERRIMAN*PBC*  *CC=  *0P=A3, JN-DBMC > 

ATTACH (MACFIL*OCMACriL, I 0*DCACSLSYS) 

ATTACH (ACSL.DCACSL*ID=DCACSLSYS) 

ACSL ( I = INPUT ) 
return*acsl.macfil. 

FTN(I«C0MPILE.R=2) 

MAP, OFF. 

ATTACH (ACSLL IB »DCACSLL IP* ID=DCACSLSYS) 
LDSET(LlR=ACSLLIB,PRESET*INOEF) 

LOAD  *LGO . 

NOGO.LGOB. 

RETURN. ACSLLIB. 

LGOR  • 

EXIT. 

0000000000000000000000 


Figure  3.  Scope  job  control  cards  for  batch  jobs 
on  the  CDC  6600. 
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MACRO  AMODPI <TN»LY*TH?) 
MACRO  REDEFINE  M,N 
LOGICAL  LY 
INTEGER  M»N»MOD» INT 
N= INT ( (THS) /3» 1416) 

M=N-N/2 

TH=THS-FLOAT (M)«6.2fl3? 

L Y =MOD ( M » 2 ) .EC.O 
MACRO  END 

Fieure  6.  l.istinR  of  AMOnPT  macro. 


MACRO  BCDCTR  (LY0.LY1  .L  Y<?  .L  Y 3 .L  YCO.LXO  *LX  1 .LX2,LX3,LXUP,LXCI  ,RST  ) 

MACRO  REDEFINE  K.KP.LHl .NABLt 
LOGICAL  LY0.LY1 .LY2,LY3.LYC0tLHl .NAOLE 

integer  k.kp.and.mino.maxo.zzl^* 

macro  RELABEL  LI  •L2.L3.L‘».LS.L8.L7,l8*L9,L10 

PROCEDURAL (L YO .L Y 1 ,L Y2 »L Y3 t LYCO=LXO ,LX 1 ,LX? ,L*3 ,LXUP .LXC I .RST ) 

IFfZZFSTIK)  .LT.0.5IGO  TO  1.2 

I F ( ZZRNFL ) GO  TO  L2 

IF  ( ZZ I CFL ) GO  TO  L6 

IF ( .N0T.NA8LEIG0  TO  L3 

IF(.NOT.TRST) (GO  TO  LI 

L8..K»MAX0(MIN0<N»7ZLSw<LX0.1.0>*ZZLSm<LX1.2.0)»ZZLS«<LX2,<i,0)»ZZLS*M... 
LX3.8.0) ) .0) 

LM 1 1 .F  ALSE . 

GO  TO  L10 
L1..K.KP-1 
IF (LXUP) K=KP* 1 
HO. .LYCO=. FALSE. 

IF  ( .NOT.  ( (K.LE.O)  .AND..NOT  .LXUDOGO  TO  L8 

IF (LHI)GO  TO  L8 

L9..LYC0=.TRUE. 

GO  TO  L7 
L8..LYO=. TRUE . 

LY1=. FALSE. 

LY2*. FALSE. 

L Y 3* .TRUE • 

GO  TO  L5 

LA..IF(.NOT. ( (K.GE.9) .AND. LXUP) (GO  TO  L7 
IF ( .NOT.LHl (GO  TO  L9 
LYO*. FALSE. 

LY1*. FALSE. 

LYZn. FALSE. 

LY3*. FALSE. 
x = 0 

GO  TO  L5 

L7..LY0  = ANI)(K,1)  .NE.O 
L Y 1 * AND IKi?) .NE.O 
LY?*AND (K  « A) .NE.O 
L Y 3* ANO ( K . 8 ) .NE.O 
L5..KP«K 
LHI«LYCO 
L3..NABLE=LXCl 
L?. .CONTINUE 

£no 

MACRO  END 


Figure  8.  Listing  of  the  BCDCTR  macro. 
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PROGRA"  RCOCTR  TEST 


-------- 

...... 

--RROV  It 

:es  enviro) 

slMF  F 

,T 

FOR  HCDCTR 

MODULE  " 

IN  I T I AL 

LOGICAL 

l XC1 

. 1 xt  P 

. ntJM 

9 

Ap 

A 1 9 9 9 

9 

A? 

9 AT 

. Lyh 

9 

HO  . 

Rl 

9 

M P 

« <3 

. lsw 

9 

LM  I OH 

LOGICAL 

PIJMP 

CONSTANT 

OUM 

= .FALSE 

..  AO  = 

. F Al  SE  . . 

Al 

■ 

.FALSF . 

9 9 9 

• 

A?  r 

.falsf  . 

• AT  = 

.TRUE.  . 

HO 

« 

.falsf. 

9 9 9 

9 

HI  = 

.FALSE . 

9 M?  = 

.FALSE.  , 

H3 

* 

.TRUE. 

CONSTANT 

TSTP 

- ■>H. 

, Dump 

= .FALSE. 

C INTERVAL 

CINT 

=:  1 .0 

AlOOPITr1'* 

I alg 

S U 

NSTEPS  NSTP  = 

1 

MAXTERVAL 

«AXT 

= 1.0 

END 

DYNAMIC 

NUMHEP  = AMAKl  (ftMlNl  (9*).  , RSW (I XO,  1 . .0.  > ♦ RSW  ( L X 1 . 7 . , o . ) 

♦ RSW (LX?, A. .0. > « »Sw(LX3.P..O.)  ... 

♦ RSW (LY0.10..0.)  . RSW (l Yl  ,?0.  .0. ) 

♦ RSW<LY?,a0.,0.)  ♦ RSW (LY3.R0. .0.) ) * 0.) 

DERIVATIVE 

L *C I = LSW(T.LT.3P.,  CLOCK (3.0) .AND. (T.OT4P.0)  . CLOC*<?.0)> 

LXUP  = LSW  ( T .C-T.  3*..  .FALSE..  .TRUE.) 

PROCEDURAL  <LYH=LHIGM> 

»CDCT»<L*0.  LXl.  LX?,  L X 3.  LXRrAO.  Al.  A?,  A9,  LXIIP,  LXCI,  LhIGh) 
RCOCTRilYC.  L Y 1 « IY?,  L Y 3 , LYP.PO , HI.  R?.  R3.  LXUP.  LXp.AND.  ... 
LXCI,  LHir.H) 

LKICH  = L xh  .tun.  < ( <L  YO.  AND.l  Y3)  .and.  LXi)P>  .OR.  .NOT.ILYO  ... 

.OP. I Yl ,OR.LY?.0R.L Y3.OR.LXUP) ) 

END  % "PROCEDURAL" 

END 

TfPMKT  ,OT.  TSTP) 

FND 

TERMlNAi 

IE (DUMP)  CALL  OFRUG 

END 

END 

ooo oooo ooo oooooooooooe 

PCEPAR  T .NUMHFR.LXCI . l x*  ip  *l yh .l  xh.lh I GH 
SET  l)U“P  = .TRUE. 

SET  TSTPsAO. ,NDRUG=S 

PROCED  GO 

START 

FRINT  "ALL" 

END? "OF  00" 

GO  « REINIT  t SFT  TSTF  = ?a(1.  % GO 
STOP 

0000 00 00000000 000000 00 
000 00000000000000000 CO 


Figure  9.  ACSL  test  program  and  associated  run 
time  commands  for  BCDCTR. 
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Figure  10. (a)  Test  results  for  BCDCTR. 
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Figure  10. (b)  Test  results  for  BCpCTR 
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Figure  10. (d)  Test  results  for  BCDCTR. 
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Figure  10.  (e)  Test  results  for  BCDCTR. 
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MACRO  81NCTR ILYO.LYl ,LY£.LY3*LYCO.LXO.LX1 .LX2,LX3.LXUP,lXCI. RST> 

MACRO  REDEFINE  K .KP ,LH 1 .NABLt 
LOGICAL  LYO.LYl .L Y2 .L Y 3 ,LYCO t LH I . NARLE 

integer  k.*p.and.mino.maxo»zzls%i 

macro  RELABEL  LI  .L2,L3.L<*.L5.L<S.L7.LB.L9.L10 

PROCEDURAL (LYO.LYl »L Y2.L  Y3  «L  YCO=LXO«LX1 .LX2.Lx3.LXUP.LXCl.RST) 

IF (ZZFST (K) .LT.O.S1GO  TO  L2 
IF (ZZRNFL)GO  TO  L2 
IF (ZZICFL)GO  TO  L6 
IF ( .NOT .NARLEIGO  TO  L3 
IF  < .NOT • (RST ) ) GO  TO  LI 

Lfc..K-MAX0 (MINO (Y. ZZLS« (LXO. I .0) ♦ZZLStMLXl .2 . 0 ) »ZZLS* <LX2 .4 . 0 > »ZZLS*(. 
LX3.8.0) ) .0) 

LHI*. FALSE . 

GO  TO  L 1 0 
L1..K-KP-1 
IF (LXUP>K=KP*1 
L10..LYCD=. FALSE. 

IF  ( .NOT.  ( (K.LE.O)  .AND.  .NOT  .LXUP)  )G0  TO  LA 

IF (LHI)GO  TO  LP 

L9..LYC0=.TRUE. 

GO  TO  L7 

lb..lyo=.toue. 

L Y 1 * • T RUE . 

L Y2* .TRUE . 

LY  3* .TRUE . 

*=15 

GO  TO  LS 

LA.. IF (.NOT. < (K.GE.15) .AND. LXUP) ) GO  TO  L7 
IF ( .NOT .LHI ) GO  TO  L9 
LYO*. FALSE. 

LY1  = . FALSE. 

LY?«. FALSE. 

LY3*. FALSE. 

* = 0 

GO  TO  LS 

L7. .LY0=ANO(K, 1 ) .NE.O 
LY1»AND(K,?> .NE.O 
LV?*AND<K.A> .NE.O 
LY3»AND(K.B)  .NE.O 
LS. .KP*K 

lbi=lyco 

L3..NABLE=LXCI 
L2.  .CONT  INIIE 

End 

macro  end 


Figure  12.  T.istfng  of  the  BTNCTR  macro. 


PROGRAM  rinctr  test 

II.  — . 

INITIAL 


-PROVIDES  ENVIRONMENT  FOR  RINCTR  MODULE 


END 


DERIVATIVE 

LPST 


LOGICAL 

l XCI  , 

LXt'P  . DUMP 

. RO  .PI  ... 

• 

a?  , 

«3  . LRST 

. 1 Sto 

CONSTANT 

HO  = 

.FALSE.  . 

HI  = .FALSE.  , 

b2  « .FALSE. 

t 

B3  = 

.TRUE.  . 

LPST*. FALSE.  . 

TSTP  * 3P. 

t 

DUMP 

= .FALSE. 

cinterval 

CINT 

* 1.0 

ALGORITHM 

IALG 

= U 

NSTEPS  NSTP  = 

1 

maxTERVAL 

MAX  T 

= 1.0 

nc 

NUMBER  = 

AMAx 1 (AMIN1  ( IS 

..  PSW (L  YO  * 1 . • 0 

»)  ♦ RSm (LY|.?,*0.)  ... 

♦ RSN(LY?»R.«o 

. ) * HS* (LY3.P..0.)).  0.) 

F NO 


END 


= lSrf<T  .EO.  31.,  .TRI'E.,  .FAISF.) 

LXCI  = LSW(T  .GT.  2.0.  .TRUE..  .FAlSE.t 

LXUP  = LS*(T  .GT.  ?0.  . . I A(  SE • . .TRUE.) 

w I NC  TP ( L YO  , LYl.  L Y2 • LY3.  L YCO  = RO.  Rl«  R?«  LXUP.  LxCI 

. LPST) 


TEPVTIT  .GT.  TSTP) 


TERMINAL 

IE  (DUMP)  CALL  DERI  iG 

END 

END 

OOOOOOOOOOOOOOCOOOOOCO 

PREPAR  T .NUMBER. LXCI  . I xi 'P  ,L  YCO  . I. pS  T 

SET  NOHUG  = 10 

PROCED  GO 

START 

PRINT  “ALL" 

ENOVOE  GO" 

GOSPF I N 1 T*SE  T TSTP=<.S.«GA 
STOP 

0000000000000000000000 
00000000000 000 000000 AO 


Figure  13.  ACSI.  tesf  program  and  associated  run  time 
commands  for  BINCTR. 
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LINf 

T 

NUMBER 

LXCI 

LXUP 

LYCO 

L»ST 

0 

0. 

8.00006 

♦ 00 

F 

T 

F 

F 

1 

1 .0900E*00 

8. OOOOE 

♦ 00 

F 

T 

F 

F 

2 

2 .00006*00 

8. OOOOE 

*00 

F 

T 

F 

F 

3 

3 . OOOOE ♦ 00 

8.00006" 

*00 

T 

T 

F 

F 

4 

4.00006*00 

8.00006 

♦ 00 

T 

T 

F 

F 

5 

S . OOOOE  *00 

9. OOOOE 

♦ 00 

T 

T 

F 

F 

6 

6 . OOOOE  » 00 

1. OOOOE 

♦ 01 

T 

T 

F 

F 

7 

7 . OOOOE ♦ 00 

1.10006 

♦ 01 

T 

T 

F 

F 

f> 

«.00OOE»00 

1.20006 

*01 

T 

T 

F 

F 

9 

9 . OOOOE  *00 

1 .30OOE 

♦ 01 

T 

T 

F 

F 

10 

1 .O000E»01 

1.4000E 

*01 

T 

T 

F 

F 

1 1 

1 . 1 OOOE  » 0 1 

1 .50006 

♦ 01 

T 

T 

T 

F 

12 

1.20006*01 

0. 

T 

T 

F 

F 

13 

1 . 70006*01 

1.00006 

♦ 00 

T 

T 

F 

F 

14 

1 .4000E*01 

2.00006 

♦ 00 

T 

T 

F 

F 

IS 

1 .5P09E»01 

3.00006 

♦ 00 

T 

T 

F 

F 

IS 

1 .6000E *01 

4.00006 

♦ 00 

T 

T 

F 

F 

17 

1.70006*01 

5.00036 

♦ 00 

T 

T 

F 

F 

IB 

l ,oroot«oi 

6.00006 

♦ 00 

T 

T 

F 

F 

19 

1 .90006*01 

7 . 00  OOF 

♦ 00 

T 

T 

F 

F 

20 

2 .OOOOE ♦ 0 I 

8.00006 

♦ 00 

T 

T 

F 

F 

21 

2 . 1 OOOE ♦ 0 1 

9.000PF 

♦ 00 

T 

F 

F 

F 

•22 

2.20006*01 

8.COOOE 

♦ 00 

T 

F 

F 

F 

23 

2.30006*01 

7.CC0OE 

• 00 

T 

F 

F 

F 

24 

2.40006*01 

6. OOOOE 

♦ 00 

T 

F 

F 

F 

2S 

2.50006*01 

5.0000C 

♦ 00 

T 

F 

F 

F 

26 

2.60006*01 

4.0000* 

• 00 

T 

F 

F 

F 

27 

2.70006*01 

3.00006 

♦ 00 

T 

F 

F 

F 

23 

2.B0006 .01 

2.00006 

♦ 00 

T 

F 

F 

F 

29 

2.90006*01 

1.00006 

♦ 00 

T 

F 

F 

F 

30 

3.00006*01 

0. 

T 

F 

T 

F 

31 

3.100C6 *01 

1.SC0C6 

♦ 01 

T 

F 

F 

T 

32 

3 . ?0GCE ♦ 0 1 

B.OOPOf 

♦ 00 

T 

F 

F 

F 

33 

3.70006*01 

7.0COOK 

♦ OJ 

T 

F 

F 

F 

34 

3.40006*01 

6.0000C 

♦ 00 

T 

F 

F 

F 

35 

3.50006*01 

5.00006 

♦ 00 

T 

F 

F 

F 

36 

3.60006*01 

4 . 0 0 0 OF 

♦ 00 

T 

F 

F 

F 

37 

3.70006  *01 

3.00006 

♦ on 

T 

F 

F 

F 

3P 

3.BOOCE.01 

2.00006 

♦ no 

T 

F 

F 

F 

39 

3.90006*01 

l .00006 

♦ 00 

T 

F 

F 

F 

LINF 

T 

NUMBER 

LXCI 

LXUP 

LYCO 

LRST 

0 

3.90006*01 

1 .OOOOE 

♦ 00 

T 

F 

F 

F 

1 

4.0000E*01 

0. 

T 

F 

T 

F 

2 

4.10006*01 

I .50006 

♦ 01 

T 

F 

F 

F 

3 

4 . 2000E ♦ 0 1 

1 .40006 

♦ 01 

T 

F 

F 

F 

4 

4.70006*01 

1.30006 

♦ 01 

T 

F 

F 

F 

5 

4.40006*01 

i .?rcoc 

♦ 01 

T 

F 

F 

F 

6 

4.5000E*01 

1 . 1 '1906 

♦ 1)1 

T 

F 

F 

F 

7 

4, 40006 *01 

l . r o o r r 

*01 

T 

F 

F 

F 

STOP 


Figure  14.  Test  results  for  BINCTR 


MACRO  CLOCK (LY » R ) 

DIFF (LY=LPULSE  <0.0*Pt(R)/10.0>) 
MACRO  END 


Figure  15.  Listing  of  CLOCK  macro. 


PROGRAM  CLKAPLSAMODASWASTP  TEST 


» PROVIDES  ENVIRONMENT  FOR  CLOCK,  LPULSE, 

••  MODE.  SWITCH.  AND  LSfEP  TESTS 

INITIAL 

LOGICAL  DUMP 

LOGICAL  LX  , Ly5 


CONSTANT  X = 1.0  . PERIOD  = 2.0  . TZ  - ?.0 

, P = 4.0  . W = 1.0 

CONSTANT  DUMP  =. FALSE..  TsTP  « 15.  . RmN  • l.E-30 

CINTERVAL  CINT  = 1.0 

ALGORITHM  I ALG  = 4 

NSTEPS  NSTP  * 1 

MAXTERVAL  maxt  = 1.0 

END 

DYNAMIC 

DERIVATIVE 

INPUT  = STEP (1.0) 

LX  * STEP (1.0)  - STEP (3.0)  .GT.  O.S 

CLOCK  <LY1=PERI00) 

LPULSE  <LY2«TZ.P.W) 

MOOE  <LY3=X) 

SWITCH<Y4=LX, INPUT) 

LYS  = (.NOT.  (LSTEPO.O)  .ANO.LSTEP  (5.0)  ) ) .OR  .LSTEP  < 7. 0) 

END 

TFRMT ( T .GT.  TSTP) 

END 

TERMINAL 

IF (DUMP) CALL  DERUG 

END 

END 

0000000000000000000000 

PREPAR  T.LY1.LY2.LV3.Y4. INPUT. I x.LYS 

SET  NDBUG  = 8 

SET  TSTP«2.0SSTARTSPRINT  "ALL"$RE INI TtSET  TSTP«P .0»ST ARTtPR 1ST  "ALL" 
STOP 

0000000000000000000000 

0000000000000000000000 


Figure  16.  ACSL  test  program  for  the  CLOCK,  LPULSK,  MODF., 
SWITCH  and  LSTEP  mcaros. 
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THS=  |*THDOTdt+THICP 


PTR(X,V=R.THSI 
AMODPKTHP.LY  = THS) 
TH  = RSW(LV.+THP.-THPI 


C 1*0 


Fipure  18.  Flow  diapram  for  CFTR. 

MACRO  CPTR  <X,Y,TH,LY.R«THOOT»THIc> 
macro  redefine  THS«THP,THICP 
MACRO  RELABEL  LI.L2 
PROCEDURAL  <THICP=THIC> 

IF ( . NOT . ZZRNFL ) GO  TO  LI 

TMICP=THS 

GO  TO  L2 

L1..IF(.N0T.Z7ICFL>G0  TO  L2 
AMOOPI (THICP.LY=THIC) 

L?.. CONTINUE 
END 

THSsINTEGITHDOT, (THICo) ) 
P7R(X.YrR,THS) 

AMOOPI (THP,LY=THS) 
TH=PSW(LY.*THP.-THP) 

MACRO  END 

Fipure  19.  Listinp  of  CPTR  macro. 
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PROGRAM  CPTR  TEST 

h PROVIDES  ENVIRONMENT  FOR  CPTR  TEST 

INITIAL 

LOGICAL  DUMP 

CONSTANT  THIC  = 0.1745,  Ol)MP  *. FALSE.,  TSTP  * 5. 

. THDOT  = 1.745,  R = 1.  . RADEO  =:  57.3 

CINTERVAL  CINT  = 0.1 
ALGORITHM  I ALG  = 4 
NSTEPS  NSTP  = 1 
MAXTERVAL  MAXT  = 0.1 

END 

DYNAMIC 

DERIVATIVE 

CPTR(X .Y.TH,LY=R. THDOT, THIC) 

LYG  = RSVMLY.  1.0.  0.0) 

END 

THG  = TH*RAOEG 

THOOTG  = THDOT*RADEG 

THICG  = THIC*RADE0 

TERMT (T  «GT . TSTP) 

END 

TERMINAL 

IF(DUMP)CALL  DERUG 

END 

END 

oooooooooooooooooooooo 


Figure  20.  ACSL  test  program  for  CPTR. 


SET  PRN«9.CMD=DIS 

SET  PRNPLT*. FALSE. .CALPLT*. TRUE. ,TTLCPL». TRUE. 
PREPAR  T.X.Y.THG.LYG 

SET  TITLE*»17  APRIL  78-CPTR  TESTER". TSTP»15. 
STARTSPLOT  "XHI"«15..X.Y 
PLOT  THG.LYG 

DISPLY  THICG.THOOTG.THG.X.Y »T»R 

SET  THDOT »-l  .745SSTARTSPL0T  "XHIM.15..X.Y 

PLOT  THG.LYG 

DISPLY  THICG. THOOTO.THG.X.Y.T.R 
PRINT  "ALL” 

STOP 


Figure  21.  ACSL  run  time  commands  on  the  TEKTRONIX  4014. 
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THG  • 10' 


2.00 


1 20 


0 40 

CM 


-0  40 


-1.20 


-2  00 


0 ANO  0 COINCIDE 


T 


Figure  25.  THC  and  1.Y0  for  THDOT  = -1.745  rad/sec. 
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ZZRNFl 

7 


T 


THICF'THS 


Figure  26.  Flow  diagram  for  CRA. 

MACRO  CRA  (X.Y,TH,i_Y.U,V.THDOT»THIC) 
WACRO  REDEFINE  THS»THP  »THICP 
MACRO  RELABEL  H»L? 

PROCEDURAL (THICP=THIC) 

IE ( .NOT . 77RNEL ) GO  TO  l 1 

THICP  = Th«; 

GO  TO  L2 

L1..IE(.nOT./ZICFL»GO  to  l? 

AMODPI (ThICP*LY=THIC) 

L?. .CONTINUE 
END 

THS=INTEG (THOOT, (THICPII 
RA (y «Y=U.V.THS) 

AMODPI < T HP  » L Yr THS) 

THxRSW (LY ♦ ♦THP  »-THP) 

MACRO  END 

Figure  27.  Listing  of  CRA  MACRO. 


i 
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PROGRAM  CRA  AND  RA  TEST 

PROVIDES  ENVIRONMENT  FOR  CRA  ♦ RA  TEST  •' 

INITIAL 

LOGICAL  DUMP 

CONSTANT  THIC  * 0.1745*  DUMP  a.FALSE. « TSTP  * 5.  ••• 

* THDOT  » 1.745*  U » 1.  . V • I. 

• RADEG  = 57.3 
CINTERVAL  CINT  » 0.1 
ALGORITHM  IALG  a 4 
NSTEPS  NSTP  * 1 
MAXTERVAL  MAXT  a 0.1 

END 

DYNAMIC 

DERIVATIVE 

CRA(X.Y.TH*LY»U.V.TND0T.THIC) 

LYG  3 RSV (LY  * l.O*  0.0) 

END 

THG  = TH*RADEG 

THICG  = THIC*RADE6 

THDOTG  a thdot«radeg 

TERMT (T  .GT.  TSTP) 

END 

TERMINAL 

IF (DUMP) CALL  DEBUG 

END 

END 

000000000000000000000000 

0000000000000000000000 

Figure  28.  ACSL  test  program  for  CRA. 

| 
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SET  PRNa9*CMD*DIS 

SET  PRNPLTa. FALSE. *CALPLT*. TRUE., TTLCPLa. TRUE. 

PREPAR  T ,X,Y»TMG,LY6 

SET  TITLE»m17  APRIL  78-CRA  ♦ RA  TESTER"*TSTP«15. 
STARTSPLOT  "XHl,,al5.  *X  . Y 
PLOT  THG, LYG 

OISPLY  THICG, THDOTG, THG*X*Y,T*U,V 

SET  THOOTa-1 .745SSTARTSPL0T  "XHI 15* , X ,Y 

PLOT  THG, LYG 

OISPLY  THICG, TH00T6, THG, X,V,T,U,V 
STOP 

Figure  29.  ACSL  run  time  commands  on  the  TEKTRONIX  4014. 
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Figure  34.  Flow  diagram  for  ('.RTF. 
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macro  crtp(R*th,thdot*ly*x.y> 

macro  redefine  ths*thi »oldly*thl •tl*thdotl*tmotic 

macro  relarel  li  .L2.l3.l4,l5*l<s.L7 

LOGICAL  LY«OLDLY 

Constant  thdtic=o.o 

RTP(R»THS=X*Y) 

procedural <th.thdot.ly=ths> 

IF(ZZFST(THS> .GE.0.5)GO  To  Ll 
IF(ABS(THS-THL)  ,GE.3.1M6)LY=.N0T.0LDLY 
GO  TO  L2 

Ll..IF(.NOT.ZZRNFL)GO  TO  L3 

Ly*OLDLY 

ThI=THL 

Go  TO  LA 

L3..IF(.N0T.ZZICFL)G0  TO  L5 
LY=.TRUE. 

GO  TO  L6 

L5..IF(ABS(THS-THL) .GE. 3. 1416)00  TO  L7 

L6..THI=THS 

GO  TO  L4 

L7..LY=.N0T.0LDLY 

THI*THS»SIGN(6.2832tTHL) 

L4.. CALL  Z7XD0T(TH00T*THL,TL*TM00TL=THI*THDTlC) 

oldly=ly 

Thl=THS 

L2..TH=THS 

IF ( .NOT ,LY)TH=-THS 

end 

macro  end 


Figure  35.  Listing  of  the  CRTP  macro. 


Flpure  36.  Illustration  of  Til,  LY  and  THS  for  a continuous 
monoton ical  1 v increasing  and  decreasing  THS  and  how  THT 
is  calculated  for  the  first  value  of  THS  past  a disconti- 


nuity in  order  to  obtain  the  correct  sipn  for  THDOT. 


PROGRAM  CRTP  TEST 


PROVIDES  ENVIRONMENT  FOR  CRTP  TFST 


IMTIAL 


LOGICAL 

DUMP 

. LSTEP 

CONSTANT 

UUMP=. FALSE. 

, TSTR  * 15. 

• 

RAOEG  = S7.3 

, LSTEP-. TRUE 

t 

K = 1.7A5 

CINTERVAL 

CINT  = 0.1 

ALGORITHM 

IALG  = A 

NSTfPS  NSTP  * 1 

maxterval 

MAXT  * 0.1 

END 

DYNAMIC 

DERIVATIVE 

thic  » -r 

w » 0.1 


THO  = RSWILSTFP.  K»STEP(0.0)«  K»SIN<W«T)> 

CPTR<x*Y.Tms,lXt1..Thd.Thi<-> 

CRTP  <R.TH,ThDOT.LY=X.Y> 

LYG  * RSW(LY.  1.0.  0.0) 


FND 


tmdo 

TMSO 

THG 

THOOTG 
Tf  RUT  (T 

END 


thd»radfg 

TMS*RA  DFG 
TH»PAp£G 
THDOT»PAOEG 
GT.  T ST  P ) 


terminal 

IF ( DUMP ) C ALL  DEBUG 

END 

FND 

00000000000000 00 000000 


Figure  37.  ACSL  test  program  for  CRTP. 
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SET  PRN*9.CHD=0IS 

SET  PRNPLTs.FALSE.*CAIPLT*.TBUE.,TTLCPL».TRUE. 

SET  T1TLE="9  AUG  78-CRTP  TESTER" 

PREPAR  T.X.Y.THOOTG.THG.R.THSG.THDG.LYG 

SET  LSTEP*.EALSE.»ThIC»-U.3I.K=1.4  31.TSTP  = 70. 

START5PL0T  "XHI"  = 70. » THG. THOOTG 
PLOT  THG.LYG 
PLOT  THDGtTHSG 
PRINT  "ALL" 

SET  K*1.745.THIC=-17.45SSTART$PL0T  "XHI "»70 . . THf: .THOOTG 
PLOT  THG.LYG 
PLOT  THDGtTHSG 

SET  LSTEP«=.TRUE..THIC=O.O.TSTP=?O.SSTARTSPLOT  "XHI"«?0.,TH«.ThDOTG 
PLOT  X.Y 
PLOT  THSG.LYG 

SET  K=~ 1 , 745SST ARTSPLOT  "XHI»=?0. .THG.THOOTG 
PLOT  X.Y 
PLOT  THSG.LYG 

SET  K=1.5708SSTART$PL0T  »XHI"=?0. .THG.THCOTG 
SET  K=- 1. 5708 1ST ART fPLOT  "XHI "=20 .. THG . THOOTG 
SET  TSTP*2.2SSTART SPRINT  "ALL" 

REINITSSET  TSTP=6.1*START$PRINT  "ALL" 

STOP 


Fipure  38.  ACSL  run  time  commands  on  the  TEKTRONIX  4014 


jOl  9J.OOHJ. 


01  • 9HJ. 


CRTP TESTER 


THDOTG  • 10 


THSG • 10 


200 


CRTP  TESTER 


Figure  52.  - THG  and  THPOTC  for  THD  = -1 . 57nR*STFr(0.0)  rad/sec. 
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MACRO  OfFF(Lr.LX) 

MACRO  RKOFFIhF  L XL 
LOGICAL  LXL.LY 

macro  Rf LANLL  L 1 «L2 
PROCFOL.c  AL  (LY  = LX  ) 

TM7ZFST  (I.XL  ) .LT.0.P)r.o  TO  Ll 

IF (77RNFL  >GO  TO  L 1 

IF ( .NOT.77ICFL )GO  To  | 7 

ly=.fal^f. 

lxl=.f alcf. 

GO  TO  tl 

LP.  .LY  = (L  X)  .AM).  .NOT. l XI 

lxl=lx 

L 1 . .COM  JMJF. 

FND 

macro  emt 


Figure  S4.  - Listing  of  PTFF  macro. 
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PROGBA"  OIFF  TEST 


INITIAL 

LOr-ICAL  Oi'HP 
LCF-1CAL  lx.  LY 
CONSTANT  CUPP  r.F&|_SF.. 
C1NTFRVAL  CINT  = 1.0 
ALC0P1TPM  I ALT.  r 6 
NST^PS  KSTP  = 1 
VAXTFPVAL  mAXT  : 1.0 

FMT 


PRPVIOES  ENVIRONMENT  FOR  OIFF  TFST 


IsTP  = IS.  , PMN  = l.E-30 


M 


dynamic 

DERIVATIVE 

l»  = STEP  < 1 . 0 ) - STPP(O.O)  ♦ STEP(F.O)  - STEP(R.O)  .CT.  O.S 

OIFF  (L  Y = STEP  ( 1 . 0 ) - STEP(^.O)  ♦ STEP<*.0>  - STFP(R.O)  .fT.  O.S) 

FAQ 

TERWT  ( T .C.T.  TSTP) 

END 

TERMINAL 

IF  (Cump  ) C all  nE°ur. 

F*  D 
ENO 

000000000000 0C0C00O00C 

PREPAR  T.LX.LY 
SFT  00^10=0 

SF  T TSTP  = S . ‘ST ART«PR!x T '•  ALL"SS E I M T SSF T TSTP= 1 0 . *ST ARTSPR I N T "ALL"SSTOP 

OOOOOOOOOCCOOOOCOFAOfr 

00000000000000000 COP 00 


Figure  55.  - ACSL  test  program  and  run  time  commando  for  DIFF. 
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line  T lx 

0 0.  F 

1 1 • 0000E ♦ 00  T 

?.  ^ .OCOvF  *00  T 

3 3. 000 OF *00  T 

A A.0000E*00  F 

S S.0000E*0O  F 

* F.0000E*00  T 


LINE  T LX 

0 *.oooor*oo  t 

1 7. 000 OE ♦00  T 

? ».  00  0 Or.  *0Q  T 

3 *5. 000  OE  ♦CO  F 

A 1.030 Of *01  F 

5 1 . 1 0 OOt ♦ 0 1 F 


r 


MACPO  DLYFF <LY. LX) 

MACRO  REDEFINE  OLDX 
MACRO  RELABEL  L1»L2.L3 
LOGICAL  OLDX.LY 
PROCEDURAL (LY=LX) 

IF (ZZFST (OLDX) *LT.0.5)G0  TO  LI 

IF (ZZRNFL)GO  TO  LI 

IF  (.NOT.ZZICFDGO  TO  L3 

ly=lx 

GO  TO  L2 

L3..LY=0LDX 

L2 . « OLDXsLX 

LI. .CONTINUE 

END 

< MACRO  END 


Figure  58.  - Listing  for  DLYFF  macro. 


program  nLYFF  TFST 

" PRO V I OES  ENVIRONMENT  FOR  DLYFF  TEST  •• 

INITIAL 

LOGICAL  DUMP 
LOGICAL  LA.  LY 

CONSTANT  DUMP  I. false..  TsTP  a is.  , Rmn  • 1 .E-30 

CINTERVAL  CINT  = 1.0 

ALGORITHM  I ALG  * A 

NSTEPS  NSTP  =>  1 

MAXTERVAL  MAXT  = 1.0 

END 


dynamic 

DFR I V AT  1 VE 

Lx  = STEP ( 1 . 0 > - STEP(A.O>  ♦ STEP(O.O)  - STEP(R.O)  .GT.  0.5 

DLYFF (L Y=L  X ) 

END 


TERMT ( T .GT.  TSTPI 

END 

TERMINAL 

i f ( dump ) call  oepug 

END 

END 

OOOOOOOOOOOOOOOOOOOOOO 
PPEPAR  T.LX.LY 
SET  ND0UG=fl 

| SET  TSTP*6 . SS TART  SPRINT  "ALL"*RE IN  I TSSET  TSTP* 1 O.S*TARTSPR  INT  "ALL"SSTOP 

OOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOO 


I 


Figure  59.  - ACSL  test  program  and  run  time  commands  for  DLYFF. 
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LINE  T 

0 0. 

1 1.0000E«00 

2 2 . 0000E  *00 

3 3 . OOOOE ♦OO 

4 4 . OOOOE *00 

5 5.0000E»00 

6 6 . OOOOE  *00 

7 7 . 0000E«00 


LINE 

0 

1 

2 

3 

4 


7 . OOOOE  *00 
ft.OOOOE^OO 

q.ooooe^oo 

I *0000E«0l 
1,1000E*01 


Figure  60.  - Test  results  for  DLYFF. 


MACRO  LMMINT (Y.YD.ICtll »LO*LL*UL) 

MACRO  RELABEL  L 1 tL2.L? 

MACRO  REDEFINE  YDOT  « I CP ♦ MODE  1 .MODE  t STATE  * I CPP  *FL AG*7ER0 

LOGICAL  MODE1.MODE 

INTEGER  FLAG 

CONSTANT  7ER0=0 • »FLAG=0 

YOOT=YD 

PROCEDURAL  (MODE *Y=YDOT»LI»LO»LL»llL»STATE*IC) 

IF (Z7FST(ICP) .LT.O.SIGO  TO  LI 
IF ( .NOT . 77PNFL ) GO  TO  L2 
FLAG=1 
ICP=Y 

ICPP=STATE 

L?. .CALL  ZZICS  (ICPP=ZERO) 

IF(.NOT.ZZICFL)GO  TO  L 3 

IF (FLAG.EC.O) ICP=AMAX1 (AMJN1 (IC*UL) *LL) 

MODE=. FALSE. 

GO  TO  LI 

L3.. MODE*. FALSE. 

IF (Y.LE.LL)MODE=-1..EC.SION(1.*YOOT) 

IF (Y.GE.UL)M0DE=l..rO.SIGN(l..YD0T) 

MODE  1= (l I) .AND.. NOT. <lD 

mode=mode.or.modei.or. r »lt> .and.. not* (Li > > 

IF (MOOE1) ICPP=STATE 

L1..Y=STATE*ICP-ICPP 

ENO 

ST ATE= INTEG <RS* (MODE »0.» YDOT) *0.) 

MACRO  ENO 


Figure  62.  - Listing  of  LMMINT  macro. 
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PROVIDES  ENVIRONMENT  EOR  LMMINT  TEST 


f 


PROGRAM  LMMINT  TEST 

INITIAL 

LOGICAL  DUMP 
LOGICAL  LI  * LT 

CONSTANT  DUMP  = . FALSE..  TsTP  = S.O 
CONSTANT  LL  =-?.  . UL  = P.  . WCPS  = 1.0 

. TWOPI  =8.283?.  A = 18. 

CINTERVAL  CINT  = 0.01 
ALGORITHM  IALG  = *. 

NSTEPS  NSTP  = 1 
MAXTERVAL  MAXT  = 0.01 
IC  =-A/(TWOPI*WCPS) 

END 


Dynamic 

DERIVATIVE 

YD  = A*SIN(TWOPI«WCPS*T> 

LIG  = STEP ( 1 . 0 > - STEP ( 3 .5 ) ♦ STEP(A.O) 

LTG  = STEP ( 1 . 0 ) - STEP (2.2)  ♦ STEP(?.9> 

LI  = LIG  .ST.  0.5 

LT  = LTG  ,GT.  0.5 

LMMINT (Y  = A»S IN ( TWOPI *NCPS«T» . I C .L I .LT .LL *UL ) 

END 


TFRMTIT  ,GT.  TSTP) 

END 


TERMINAL 

IE (DUMP) CALL  OERUG 

END 

END 

0000000000000000000000 

0000000000000000000000 


Figure  63.  - ACSL  test  program  for  LHHINT. 


SET  PRN=9.CML)=DIS 

SET  PRNPLT=.E..CALPLT=.T..TTLCPL=.T, 

Prepar  t.y.yu.lig.ltG 

SET  TITLE="LMMINI  TESTER-23  JAN  T9" 

Set  TSTP=0.99SSTARTSSAVE  "BASE»*REINITSSET  TSTP*1.6.NRWITG«.T. 

ACTION  "VAh"=2.191."VAL"=8.»L0Cm=NDBUG*START*REINITSSET  TSTP-2.9SSTART 

MEIN1TSSET  TSTP=3.6SACI ION  »CLEAR“SSTARTS«EINIT$SET  ) STP«5.0*START 

Set  NRWITG=.E.»PL0T  "Xrtr'sS.U.Y.YU 

^LOT  LIG 

(‘lot  LTG 

PRINT  ••ALL” 

Set  tstp=s.i.dump*.t.sstart$print  «allm 

RESTOR  ••BASE"SSET  A=6.2832»TSTP=S.O ,OUMP*.E.SSTARTJPLOT  »XMI*'»5.0.Y.YD 

print  "all" 

Stop 


Figure  64.  - ACSL,  run  time  commando  for  LMMINT  tests  on 
TEKTRONIX  4014. 

I 

90 


AD-A070  781 


UNCLASSIFIED 

2 of  2 I 

AD 

AO 70 781 


ARMY  MISSILE  RESEARCH  AND  DEVELOPMENT  COMMANO  REDSTO— ETC  F/G  9/2 
ACSL  MACROS  FOR  ANALOG  OPERATORS. (U) 

OCT  78  D B MERRIMAN 

DRDMI-T-79-7  Ml 


MACRO  LPULSE (LY.TZ*P.W) 
LOGICAL  LY 

LV=PULSE(TZ.P*W> .GT.0.5 
MACRO  END 


Figure  69.  - Listing  of  LPULSF  macro. 


MACRO  LSTEP (LY.TZ) 
LOGICAL  LY 
LY=STEP(TZ) .GT.0.5 
MACRO  END 


Figure  70.  - Listing  of  LSTEP  macro. 


MACRO  MOOE(Y.X) 
LOGICAL  Y 
Y*. FALSE. 

MACRO  ENO 


Figure  71.  - Listing  of  the  MODE  macro. 


MACRO  MONO(LY.R.LX) 

MACRO  RELABEL  L1*L2*L3*LA«L5»L6 
MACRO  REOEFINE  OLDT. OLDL Y*LXP 
LOGICAL  OLDLY *LY »LXP 
PROCEDURAL (LY*R»LX) 

IF  (ZZFST  (OLDLY)  .LT.f).*5)GO  TO  LI 
IF(ZZPNFL)GO  TO  LI 
IF(.NOT.ZZICFL)GO  TO  L3 
ly=lx 

IF (.NOT.LYIGO  TO  L2 
OLDT=T 
GO  TO  L2 
L3..LXP*LX 

IF(.NOT.OLDLY)GO  TO  LA 
IF <T-OLnT.LT. (R) )GO  TO  L5 
LA..IF(LXP)GO  TO  L6 
LY=. FALSE. 

GO  TO  L2 
L6..01DT=T 
L5 . .L Y= . TRUE • 

12. . OLDLY=LY 
LI. .CONTINUE 
FND 

MACRO  END 


Figure  73.  - Listing  of  the  MONO  macro. 


PROGRAM  MONO  TEST 

PROVIOES  ENVIRONMENT  FOR  MONO  TFST 

INITIAL 

LOGICAL  DUMP 

LOGICAL  LX  * LY 


CONSTANT  DUMP  =. FALSE..  TSTP  * 15.  . RMN  ■ l.E-30 

CONSTANT  R =2. 

CINTERVAL  CINT  = 1.0 
ALGORITHM  IALG  = A 
NSTEPS  NSTP  = 1 
MAXTERVAL  MAX T = 1.0 

END 


DYNAMIC 

DERIVATIVE 

LX  = STEP ( 1 .0 ) - STEP (2 . 0 ) ♦ STEP(G.O)  - STEP<8.0>  ♦ 

STEP  ( R . 0 ) .C-T.  O.S 
MONO (LY=R.LX> 

END 


TFRMT ( T .GT.  TSTP) 

END 

TERMINAL 

IF (DUMP) CALL  DEBUG 

ENO 

END 

0000 000 OOC 00 0000000000 
PREPAR  T.  lx.  LY 

SET  TSTP= 7.0SSTARTSPRINT  ••  ALL"SRF I N I TSSFT  TSTP* 1 2 . »ST ART$PR INT  "ALL" 
STOP 

0000000000000000000000 

0000000000000000000000 


Figure  74.  - ACSL  teat  program  and  run  time  commands  for  MONO. 


98 


T 


LX 


LIME 

0 0.  F 

1 1 «0000E*00  T 

2 2.0000E»00  F 

3 3.0000E«00  F 

4 4.0000E*00  F 

5 S.0000E«00  F 

6 6.0000E»00  T 

7 7.0000E«00  T 

8 8.0000E«00  F 


LINE  T LX 

0 8.0000E*00  F 

1 R.OOOOE^OO  T 

2 1 .OOOOE*01  T 

3 1.1 0OOE»Ol  T 

A 1.2000E»01  T 

5 1 . 3000E*01  T 


Figure  75.  - Test  results  for  MONO 


MACRO  RA (X.YtlUVtTH) 
MACRO  REDEFINE  C*S 
S*SIN(TH> 

C=COS(TH) 

X=U*C-V«S 
Y=U*S*V*C 
MACRO  END 

| 


Figure  76.  - Listing  of  the  RA  marro. 


MACRO  SHIFT <LYO.LYl,LY2.LY3,NICtLXSH.LXSI) 

MACRO  RELABEL  L1*L2.L3 

MACRO  REDEFINE  LLYO.LLY1 *LLY2*LLY3tK 

LOGICAL  LY0tLYl.LY2.LY3 

LOGICAL  LLYO.LLY1 .LLY2.LLY3 

INTEGER  K ♦ AND 

PROCEDURAL <LY0.LYl.LY2.LY3=NIC.LxSH.LXSI> 
IF(ZZFST(K) ,LT.0.5)G0  TO  Ll 
IF(ZZRNFL)GO  TO  Ll 
IF(.NOT,ZZICFL)GC  TO  L2 
K=NIC 

LYO=AND (K  « 1 ) .NE.O 
LY 1*AND (K  »2) .NE.O 
LY2*AND(K.4) .NE.O 
LY3=AND (K.8) .NE.O 
GO  TO  L3 

L2..IF (.NOT. (LXSH) )GO  TO  Ll 
LY0*LXSI 

lyi*llyo 

LY2=LLY1 
LY3=LLY2 
L3..LLY0=LY0 
LLY1*LY1 
LL Y2=L Y? 

LLY3=LY3 
Ll.. CONTINUE 
END 

MACRO  END 


Figure  78.  Listing  of  the  SHIFT  macro. 


PROGRAM  SHIFT  TEST 

PROVIDES  ENVIRONMENT  FOR  SHIFT  TEST 

INITIAL 

LOGICAL  DUMP 

LOGICAL  LASH  * LXSI  , LYO  . LY1  . LY? 

. LY3 

INTEGER  NIC 

CONSTANT  DUMP  =. FALSE..  TsTP  a 15.  , RMN  » l.E-30 

CONSTANT  NIC  = 15 

CINTERVAL  CINT  = 1.0 

ALGORITHM  I ALG  = 4 

NSTEPS  NSTP  * 1 

MAXTERVAL  MAXT  = 1.0 

END 


k 


DYNAMIC 

DERIVATIVE 

LXSI  * STEP (2.0)  - STEP ( 3 .0 ) ♦ STEP(4.0)  - STEP(5.0)  ♦ 
STEP (9.0)  .GT.  0.5 

LXSH  = STEP ( I . 0 ) - STEP (6.0>  . STEP(7.0>  *GT.  0.5 
SH IFT (LYO. LY1,LY2.LY3=NIC,LXSH, LXSI) 

END 

TERMTIT  .GT.  TSTP) 

END 


terminal 

IF (DUMP) CALL  DEPUG 

ENO 

END 

0000000000000000000000 

PREPAR  T.LXSI .LXSH.LY3.LY2.LY1 .LYO 

SET  TSTPaS. OSS TART  SPR I NT  "ALL"SRF IN  I TSSET  TSTP*1 1 .OSSTARTSPRINT  "ALL” 
STOP 

0000000000000000000000 

0000000000000000000000 
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MACRO  SPTEFF (LY.LIC,LS»LR*LE) 

MACRO  RELABEL  L1.L2.L3 
MACRO  REDEFINE  OLDY.LSP.LRP 
LOGICAL  OLDY,LY,LSP,LRP 
PROCEDURAL <LY=LIC.LS«LB»LE> 

IF (72FST (OLDY) ,LT.0.5)GO  TO  LI 
IF (ZZRNFL)GO  TO  LI 
IF  (.N0T.7ZICFDG0  TO  L3 
LY=LIC 
GO  TO  L2 

L3..IFI.N0T.  (LE)  ) GO  TO  L 1 

LSP=LS 

LPP=LR 

IF (LSP ) L Y= .TRUE . 

IF (LRP)LY=. FALSE. 

IF ( ( .NOT. LSP) .AMD. (.NOT.LRP) )LY=OLDY 

IF (LSP.AND.LRP)LY=.NOT.OLOY 

L2..0L0Y=LY 

LI. .CONTINUE 

END 

MACRO  END 


Figure  82.  Lifting  of  the  SRTF.FF  macro. 
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PROGRAM  SRTEFF  TEST 


PROVIDES  ENVIRONMENT  FOR  SRTEFF  TFST 

INITIAL 

LOGICAL  DUMP 

LOGICAL  LY  * LIC  . LS  . LR  . IE 


CONSTANT  DUMP  =. FALSE..  TSTP  ■ IS.  ♦ RMN  * 1 .E-30 
CONSTANT  LIC  = .FALSE. 

CINTERVAL  CINT  = 1.0 
ALGORITHM  I ALG  = U 
NSTEPS  NSTP  = 1 
MAXTERVAL  MAXT  = 1.0 

END 

DYNAMIC 

DERIVATIVE 


LE 

= STEP (1.0)  - STFP(fi.O) 

♦ STEP ( 1 0 . 0 ) 

.GT.  0.5 

LP 

= STEP  (0.0)  - STEP  (*>.0) 

» STEP ( 7 .0 ) 

- STEP  < 1 ?»0 ) . 

LS 

= STEP (0.0)  - STEP ( 3 .0 ) 

♦ STEP (6.0) 

- STEP (11.0)  ♦ 

STEP (13.0)  .GT.  0.5 

SRTEFF  <LY=LIC.LS.LR.LE> 


END 

TERMTIT  .GT.  TSTP) 

ENO 

TERMINAL 

IE (DUMP) CALL  OEPUG 

END 

END 

0000000000000000000000 
PREPAR  T.LE.LRtLS.LY 

SET  TSTPaR. OSS TART SPRINT  "ALL"SRf INITSSET  TSTP«1 A.OSSTARTSPRINT  "ALL 
STOP 

SRTFF(LY«LIC.LS.LR) 

0000000000000000000000 

0000000000000000000000 


Figure  83 


ACSL  test  program  and  associated  run  time  commands 
for  SRTEFF  and  SRTFF. 
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Figure  8*u  Test  results  for  SRTEFF. 
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MACRO  SRTFF(LY«LIC»LS.LR) 

MACRO  RELABEL  L1«L2,L1 
MACRO  REDEFINE  OLDY,LSP.LRP 
LOGICAL  OLDY*LY»LSP»LRP 
PROCEDURAL (LY=LIC.LS«LP) 

IF (7ZFST (OLOY) ,LT.0.5)G0  TO  LI 
IF (7ZRNFL)GO  TO  LI 
IF  (.N0T.Z7ICFDG0  TO  Ll 

ly=lic 

GO  TO  L? 

L3..LSP=LS 

LRP=LR 

IF (LSP) L Y= . TRUE . 

IF (LRP)LY=. FALSE. 

IF (( .NOT. LSP) .ANO. < .NOT.LRP) >LY=OLOY 

IF (LSP.AND.LRP)L Y=. NOT. OLOY 

L?..OLDY=LY 

Ll.. CONTINUE 

ENO 

MACOO  END 
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Figure  86.  - Test  results  for  SRTFF. 


MACRO  SWITCH (Y *LX*X) 
Y=RSW(LXtX,0.0) 

MACRO  ENO 


Figure  87.  - Listing  of  the  SWITCH  macro 


110 


■ 


MACRO  TIMER(LY*LYCO»NIC*LXCI »LS*LR) 

MACRO  RELABEL  L1*L2»L3*L4 

MACRO  REDEFINE  KPtSTAPTD.K 

INTEGER  KP.K 

LOGICAL  LY*LYCO.STAPTO 

PROCEDURAL  <LY*LYCO»NlC  «LXCI*LS»LR) 

IF (ZZFST (KP) .LT.0.5)G0  TO  LI 
IF (ZZRNFL)GO  TO  LI 
LYCO». FALSE. 

IF(.NOT.(ZZICFL.OR.(LR).OR..NOT.((LS>.OR.STARTD)))00  TO  L2 
LY«. FALSE. 

L3..START0a. FALSE. 

kp«nic 

GO  TO  LI 

L2..STARTD».TRUE. 

LY*.TRUE. 

IF(.NOT.(LXCI))GO  TO  Ll 
K*KP- I 

IF(K.EO.O)GO  TO  L4 

kp=k 

GO  TO  Ll 
L4..LYCO*.TRUE. 

GO  TO  L3 
Ll.. CONTINUE 
END 

MACRO  ENO 


I 

Figure  89.  - Listing  of  TTMKR  macro. 
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PROGRAM  TIMER  TEST 

" PROVIDES  ENVIRONMENT  FOR  TIMER  TEST 

INITIAL 

LOGICAL  DUMP 

LOGICAL  L Y • LYCO  , LXCI  . LS  . LR 

INTEGER  MC 


CONSTANT  NIC  * 7 

CONSTANT  DUMP  a .FALSE • * TsTP  ■ 15.  . RMN  • l .E-30 

CINTERVAL  CINT  « 1.0 
ALGORITHM  IAL6  a <• 

NSTEPS  NSTP  * 1 
MAXTERVAL  MAX T a 1.9 

END 

DYNAMIC 

DERIVATIVE 

LXCI  a STEP (1.01  - STEP (5.0 ) ♦ STEP<7.0»  .GT.  0.5 

LS  * STEP (1.0)  - STEP (2.0)  ♦ STEP<1?.0)  .GT.  0.5 

LR  a STEPU5.0)  - STEP  (16.0)  .GT.  0.5 

TIMER(LY.LYCO=NIC.LXCI.LS,lR> 

END 

TERMTU  .GT.  TSTP) 

END 

TERMINAL 

IF (DUMP) CALL  0ERU6 

END 

END 

0000000000000000000000 
PREPAR  T.LY. LYCO. LXCI. LS.LR 

SET  TSTP= 1 e.SST ARTSPR INT  " ALL"SPE INI TTSET  TSTP»23.9.NO0UGa5«STAPT 
PRINT  "ALL" 

STOP 

0000000000000000000000 

0000000000000000000000 


Figure  90.  - ACSL  teat  program  and  run  time  commands  for  TIMER. 
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Figure  91.  Test  results  for  TIMER. 
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